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The infra-red absorption spectrum from 2 to 14.54 was determined for cyclohexane vapor. 
The number of bands found in the 700-1350 cm™ region confirms the chair structure (sym- 
metry D32) as the most stable at room temperature. The band envelopes are discussed and 
compare well with the theory of Gerhard and Dennison. 





I. INTRODUCTION 


HE infra-red absorption spectrum of cyclo- 
hexane vapor has an important bearing on 

the question of the configuration of the cyclo- 
hexane molecule. The evidence in favor of the 
“chair’’ form (symmetry D3), and hence of the 
staggered configuration of groups about a C—C 
single bond, has been summarized by Schomaker 
and Stevenson! and by Kohlrausch and Wittek.? 
Langseth and Bak* pointed out that if the 
opposed configuration of groups about a C—C 
single bond were the more stable, the cyclohexane 
molecule would assume the “‘cradle’”’ form (sym- 
metry C2,) or the planar form (symmetry Dg,) 
or some compromise between the two. The exact 
form in this case would depend on the relative 
magnitudes of the potential restricting rotation 
about the C—C bond and the energy of distor- 
tion of the tetrahedral valence angles of the 


carbon atoms. The latter authors also presented’ 


data drawn from the Raman spectrum of liquid 
cyclohexane to indicate that the structure is the 
planar one. Kohlrausch and Wittek? called 


'V. Schomaker and D. P. Stevenson, J. Chem. Phys. 8, 
637 (1940). 

2K. W. F. Kohlrausch and H. Wittek, Zeits. f. physik. 
Chemie B48, 177 (1941). 

3 A. Langseth and B. Bak, J. Chem. Phys. 8, 403 (1940). 


attention to an error in the polarization data 
used by these authors and showed that the 
correct polarization data lead unambiguously to 
the chair structure, assuming the usual selection 
rules to hold in the liquid. As shown below, this 
is also the conclusion drawn from the infra-red 
spectrum of the vapor, where the selection rules 
hold exactly. 

The Raman spectrum of liquid cyclohexane is 
recapitulated in Table I. The values shown are 
averages of those given by various workers.’ 
The infra-red absorption spectrum has been 
given previously by Kettering and Sleator® but 
without any tabulated wave-length data or any 
discussion. 


II. EXPERIMENTAL DETAILS AND RESULTS 


The spectrometer employed was an automatic- 
recording rocksalt prism instrument,® capable of 


4A. Petrikaln and J. Hochberg, Zeits. f. physik. Chemie 
B3, 217 (1929); E. Canals, M. Godchot, and G. Cauquil, 
Comptes rendus 194, 1574 (1932); K. W. F. Kohlrausch 
and W. Stockmair, Zeits. f. physik. Chemie B31, 382 
(1936); R. Ananthakrishnan, Proc. Ind. Acad. Sci. 3, 52 
(1936); K. Matsuno and K. Han, Bull. Chem. Soc. Japan 
11, 321 (1936); A. I. Sidorova, Acta Phys. (U.S.S.R.) 7, 193 
(1937). 

5 C. F. Kettering and W. W. Sleator, Physics 4, 39 (1933). 

6 To be described in a forthcoming article in The Review 
of Scientific Instruments. 
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resolving six reciprocal centimeters at eleven mu. 
The total absorption path used was 58 cm. The 
spectrum was obtained at a pressure of 40 mm 
Hg, and at regions of heavy absorption at 5 mm 
Hg also. 

The cyclohexane used was a pure sample 
prepared under the direction of Professor C. E. 
Boord of the Department of Chemistry of the 
Ohio State University as part of the American 
Petroleum Institute Hydrocarbon Research Proj- 
ect in the Industrial Research Foundation of 
the University. 

The spectrum obtained (Fig. 1) is in good 
agreement with that of Kettering and Sleator, 


TABLE I. Raman and infra-red spectra of cyclohexane 
(frequencies in cm™). 








Raman (liquid) Infra-red (vapor) 





Fundamentals 382 ( 2)P 673.4 ( 2) 
426( 3)D 
802 (10) P 864.3 ( 6) 
1028( 7) D 902.9 ( 5) 
1155 ( 4) 4P 1015: ( 1) 
1266 ( 7) D 1042 ( 3) 
1261 ( 5) 
1346 ( 2) (D?; harmonic?) 
1443 ( 7) D 1456 (10) 
2853 (11) P 2865 ( 1) 
2871 ( 2) 2920 ( 1) 
2888 ( 3) D 2960 ( 9) 
2897 ( 3) resonance 
2909 ( 1) splitting 
2923 ( 8) D 
2937 ( 8) P 
Higher harmonics 762 ( 0) 1348 (0) 
1424 ( 3) 1684? ( 0) 
1464 ( 1) 1792 (0) 
2349 ( 0) 1842? (00) 
2466 ( 1) 1887 (0) 
2536 ( 0) 2004 ( 0) 
2554 ( 0) 2079 (0) 
2604 ( 0) 2179? (00) 
2629 ( 1) yi Me eb 
2663 ( 3) P 2294? ( 0) 
2698 ( 2) 2427 (1) 
2530 (0) 
2665 ( 0) 
2705 ( 5) 
3210 ( 1) 
3360 (0) 
3495 (0) 
3815 (0) 
4000 ( 0) 
4150 (0) 
4290 
te} ( 2) 








Figures in parentheses give relative intensities (based on 40-mm 
curve for infra-red). P=polarized, D=depolarized. Infra-red fre- 
quencies for \<4y rounded off to closest 5 cm™}, (?) indicates doubtful 
bands, possibly due to wiggles in energy curve. 
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but shows a little better resolution. The wave- 
lengths (Table I) are accurate to +.01y for 
A\>5yu, as indicated by a test run on ethylene 
vapor and comparison with grating data. At 
\<5y, the accuracy is somewhat lower, becoming 
about +.02u at 3u. 

The bands with centers at 7.93u, 9.60u, and 
11.57 are clearly resolved into P, Q, and R 
branches. The peak at 9.85u appears to be the 
Q branch. of a band whose P and R branches 
are obscured by the P branch of the 9.60u band. 
The band at 11y seems to be different in appear- 
ance from the other long wave bands, having 
two equally intense peaks at 10.97 and 11.075u. 
The absorption near 14.6 is seen, fram reference 
to the curve of Kettering and Sleator, to be the 
R branch of a band with Q branch at 14.85y. 


III. DISCUSSION 


The forty-eight normal modes of vibration of 
cyclohexane can be conveniently divided into 
the following sets: twelve C—H_ stretching 
vibrations, six H—C—H _ bending vibrations, 
six C—C stretching vibrations, six C-—C—C 
bending vibrations, and eighteen vibrations in 
which CH: groups twist or rock as units. This 
division, which is of course only a first approxi- 
mation to the true motions, is sufficiently accu- 
rate to enable the frequency range to be stated 
in which each set of vibrations will be found. 
For a given symmetry of the whole molecule the 
symmetry classes to which the vibrations in each 
set belong are readily found by application of 
the usual group theory methods to the set of 
internal coordinates which describe the vibra- 
tions of the required type. In Table II are given 
the frequency ranges for each set of vibrations 
and the symmetry classes into which these 
vibrations fall, for the chair, cradle, and planar 
forms. Infra-red and Raman selection rules are 
shown for each case. The frequency ranges are 
known through reasonable extrapolations from 
corresponding vibrations in the simpler molecules 
methane, ethane, and ethylene. 

The region 1350 to 700 cm— (7.4 to 14.3,), 
which includes the C—C stretching frequencies 
and CH, rocking and twisting frequencies is the 
most fruitful one to examine, since here the 
frequencies are by their nature most spread 
apart, and resolution is highest in the prism- 
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dispersed infra-red spectrum. In the observed 
spectrum there are five bands here, exactly the 
number to be predicted on the chair structure. 
The cradle and planar structures would require 
eighteen and three, respectively. Any structure 
intermediate between the cradle and planar 
would require even more, being of lower sym- 
metry than either. Hence, the infra-red spectrum 
directly confirms the chair form as the stablest 
at room temperature. Reference to the Raman 
data and to Table II shows that both in polar- 
ization requirements (the polarized frequency 
shifts 382, 802, 1155 cm being totally sym- 
metric) and in distribution of frequencies the 
chair form is indicated. 

Only one of the two expected infra-red bands 
is found in the 1400 cm™ region. However, this 
band is so strong and resolution is so relatively 
poor in this region that this may very well be a 
case of a strong band completely hiding a 


weaker one. One of the low frequency infra-red 
bands is still unobserved, lying at a wave-length 
greater than 16u. In the Raman spectrum 1400 
cm! region, only one strong line has been found, 
the depolarized 1443 cm frequency. There 
should be a polarized shift in this region, but it 
is apparently rather weak. Possibly the weak 
1424 cm line (listed as a harmonic) is this 
fundamental, although its polarization has not 
been measured because of its low intensity. The 
usual rule that the totally symmetrical (polar- 
ized) Raman shifts are more intense than the 
non-totally symmetrical seems to be violated 
here. In the 700-1350 cm region, two polarized 
and two depolarized lines are found; two more 
depolarized ones should be present but are 
apparently too weak to be observed. The shift 
at 1346 cm™ is possibly one of these. 

Any more specific assignment of the observed 
frequencies, other than the symmetry of the 


TABLE II. Symmetry classes of vibrations of cyclohexane for chair, cradle, and planar forms. 











Chair, Daa Cradle, Coy Planar, De, 
Type and range Aig Aiu Az Aw Eg Eu} Ai Az Bi Bs |Atg Atu Az Aru Big Biu Boy Bou Etg Etu Ex Eou 
12 C—H stretching 2850-3000 cm= 2 6 @32 2 32 4 2 4 2 1 Ss & F Ss st 8 1 1 1 
6 H —C—H bending 1350-1500 cm~! 1 0 O11 : Ss 2 . 2 1 1 0 Oo 90 te @ & 1 1 0 
18 CHe rocking and twisting 700-1350 cm! 1 a. 2 @ es *g * $6 § » 2 @ fs eh 2 1 1 2 
6 C—C stretching 700-1350 cm= 1 1 00 s @ 2 » ¥ 2 1 0 Oo 0 » & et € 1 1 0 
6 C =-C —C bending <700 cm- 1 0 O01 1 1 2 2 2 0 0 0 oO O 0 11 0 «0 0 1 1 
Activity: IR =infra-red, R = Raman, mstete@m @min &@© 8. Mim €§ ££ FE TCR Rm RE 
I =inactive IR IR IR 
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Raman frequencies as shown by polarization 
data, is impossible at this time. The infra-red 
band envelopes seem to be incapable of yielding 
knowledge as to the symmetry characters of the 
infra-red active frequencies, as discussed below. 
Lack of knowledge of the five inactive frequencies 
and of one infra-red and several Raman fre- 
quencies makes it of little use to attempt to 
interpret the overtones and combinations. These 
specific details of the vibrational spectrum must 
await an adequate force constant treatment for 
their elucidation. 

Using accepted carbon-carbon and carbon- 
hydrogen bond lengths (1.54 and 1.09A, respec- 
tively) and tetrahedral angles, the principal 
moments of inertia of the chair structure are 
calculated to be: 


A=B=116.8 at. wt. units- A’, 
C = 201.8 at. wt. units- A’. 


The value of the parameter B=(A/C)—1, used 
by Gerhard and Dennison’ to describe the 
asymmetry of a symmetrical top molecule, is 
— 0.421; i.e., the cyclohexane molecule approxi- 
mates somewhat a flat disk (@= —4). Reference 
to the band envelope curves for perpendicular 
bands given by Gerhard and Dennison shows 
that for this value of 8, these bands still have a 
good deal of the character of bands of a spheri- 
cally symmetrical rotator, the zero branches 
(AJ=0, AK=+1) forming a strong central 
maximum, giving the whole band the appearance 
of a parallel band. This explains the parallel 
type envelopes found in four of the five bands in 


the 700-1350 cm region of the infra-red spec-., 


7S. L. Gerhard and D. M. Dennison, Phys. Rev. 43, 197 
(1933). 
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trum, where there are to be expected four 
perpendicular and one parallel band. The ap- 
parently anomalous envelope of the 11y band 
can be interpreted as being due to a vibration 
involving a relatively large change in moments 
of inertia between the ground and first excited 
vibrational states. The resulting degrading (to 
the red, in this case) causes the Q branch at 
902.9 cm-! (11.075) to be considerably broad- 
ened; the R branch at 912.4 cm (10.97) to be 
shifted to a lower frequency than normal and to 
pile up until its maximum intensity is comparable 
to the Q branch; and the P branch to be shifted 
to lower frequency than normal and decreased 
in intensity so that it is hidden by the R branch 
of the 864.3 cm~ (11.57) band. 

Using the “‘separation’”’ function of Gerhard 
and Dennison,’ the calculated P-R separation 
of the parallel bands is 23 cm~. The perpendic- 
ular band P-R separation should be slightly less 
than this. The observed values for the 1042 cm—! 
and 864.3 cm bands (the only two which were 
free enough from interference and sufficiently 
resolved to measure accurately) are 25 cm for 
each, in very good agreement considering the 
limitations of the simple theory applied to 
actual cases. 

Recently data have been published on the heat 
capacity of cyclohexane vapor over a temperature 
range from 370° to 410°K.® Because of lack of 
knowledge of the several inactive or unobserved 
frequencies and of the degeneracies of known 
frequencies it seems improbable that any signifi- 
cant result would be obtained by a statistical 
calculation and comparison with experiment. 


8 J. B. Montgomery and T. DeVries, J. Am. Chem. Soc. 
64, 2375 (1942). 
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Intensities and polarization properties of vibrations in 
which mainly chlorine atoms are participating, are investi- 
gated theoretically for the dichlorobenzenes. It is assumed 
that the change of polarizability in these vibrations is 
connected primarily with displacements of the chlorine 
atoms and that contributions from the individual atoms 
are additive. The Raman tensors belonging to the chlorine 
vibrations are then obtained by either adding or sub- 
tracting the single C—ClI contributions with reference to 


a chosen system of axes rigidly connected with the mole- 
cule. Intensity ratios were obtained for the symmetrical 
and unsymmetrical valence vibrations and deformation 
vibrations, respectively. Very small depolarization factors 
were found for lines representing symmetrical valence 
vibrations. For deformation vibrations the simplified 
theory gives complete depolarization for all allowed 
Raman lines. The results are discussed in connection with 
existing experimental data. 





INTRODUCTION 


AMAN spectra of substituted benzenes have 
been investigated by several authors and 
correlations of the observed lines to particular 
modes of vibration have been proposed in many 
cases. In spite of the great number of vibrations 
it is possible to follow such a correlation through 
whole series of substitutions for certain vibra- 
tions. However, in general, it is extremely diffi- 
cult or impossible to find an unambiguous as- 
signment for all observed lines. Such difficulties 
were encountered in the course of an investiga- 
tion of the Raman spectra of the dichloroben- 
zenes and their depolarization factors.) ? 

In this paper, we are interested in vibrations 
of the dichlorobenzenes in which mainly chlorine 
atoms are involved and in the correlation of these 
vibrations to observed Raman lines. To supple- 
ment the existing data we have studied from a 
theoretical point of view the intensity and polari- 
zation properties of these vibrations. We have 
used the following simplified assumptions: a. 
The change of polarizability in these particular 
vibrations is connected mainly with the displace- 
ments of the chlorine atoms. b. The contributions 
to the change of polarizability from the indi- 
vidual atoms are treated as additive. 

Among the vibrations of the dichlorobenzenes 
there must occur two valence and two deforma- 
tion vibrations of the C—Cl bonds in the mo- 


1H. Sponer and J. S. Kirby-Smith, J. Chem. Phys. 9, 667 
(1941). See this paper for previous references. 

20. Paulsen, Monatsh. d. Chemie 72, 244 (1939); 
K. W. F. Kohlrausch and O. Paulsen, ibid. 268 (1939) 
-dichlorobenzene. 


lecular plane and two deformation vibrations 
perpendicular to the plane whereas there is only 
one of each kind in monochlorobenzene. The val- 
ence vibrations will have larger frequencies than 
the deformation vibrations. All of them have 
frequency values below 1000 cm. 


GENERAL THEORETICAL OUTLINE 


The induced electric moment m is given by 
m=ak, 


where E is the vector of the electric field strength 
of the incident radiation and a is the polariza- 
bility tensor. a depends on the coordinates de- 
scribing the motion of the vibrating or rotating 
molecules. We are not concerned here with the 
effect of rotation, but rather with the depend- 
ence of the polarizability on molecular vibration. 
Since the displacements of the atoms during the 
vibration are small compared to nuclear dimen- 
sions, a can be developed in powers of the normal 
coordinates g; which themselves are periodical 
with time (q;=a; cos 2mv;t, where a; is the ampli- 
tude and »; the frequency of the i-th vibration). 
Including only linear terms, we have for a typ- 
ical element a, of the polarizability tensor 


T 
Oazy 


den 


total 


T 
Ary =Aryt>, 
i 





T iR 
=A2y+)>, azy Cos 2rv;it (1) 
where 
T 
OOry iR 
a; = Ary . 





0g: 
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Here the first term ie gives the Tyndall or 
Rayleigh scattering and in the second term the 
i-th member of the sum represents the Raman 
scattering due to the i-th vibration. The contri- 
butions from the different vibrations are inde- 
pendent of each other (incoherent) and we will 
drop the index 7 when discussing the Raman 


tensor dis due to one particular vibration. 

Such a tensor (assumed to be symmetrical) 
may be split into two parts, a spherically sym- 
metrical and an anisotropic one (unless specified 


®) 

















Pag m 


Fic. 1. Modes of chlorine vibrations and symmetry 
symbols for the three dichlorobenzenes. (a) Valence vibra- 
tions; (b) Planar deformation vibrations; (c) Non-planar 
deformation vibrations. 
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, R 
otherwise, a;, means the Raman tensor a,, from 
here on) 


= ty tm =—8 +( =i) (2) 
Azy = Ary we zy Ary V3 zy f° 


Here s=(1V3)(az:+a,y+azz) is the trace of the 
tensor multiplied by (1/v3) and 6,,=1 for x=y 
and 6,,=0 for x#y. Since s is an invariant of the 
tensor, this splitting is valid for any coordinate 
system, also for the special one of the principal 
axes in which the tensor a, is diagonal. 

Intensity and polarization of the Raman light 
can be expressed with the help of these tensors. 
The intensity of the total Raman scattering in 
all directions due to one particular vibration is 
proportional to the square of the corresponding 
tensor a@;,, i.e., to the sum of the squares of all 
the components (the square is also an invariant 
of the tensor). 


2 2 2 2 2 
I~> Ory = Age t+ Ayyt+AeztAryt pi = 5?-+5?, 
zy 


where s? and 6? are the squares of the tensors 


8 b ° ° 
Qz, and a;,. However, the intensity of a Raman 
line is generaliy observed at right angles to the 
incident beam. This intensity is given by the 


expression 
Tos = 10s?+ 130. (3) 


The depolarization factor (for natural light) is 
given by the expression* 


6b? 


on 4 
pn 1052+ 136? (4) 


THEORETICAL TREATMENT OF CHLORINE 
VIBRATIONS 


With this general outline in mind we shall now 
attack the problem of the chlorine vibrations of 
the dichlorobenzenes. We shall try to obtain 
values for the depolarization factors of their Ra- 
man lines and also an estimate of the intensity 
ratios between the lines belonging to the two 
valence vibrations in the molecular plane, to the 

* The formula frequently found in the literature 

6b” 
Pn 45a2+ 70" 
follows from the one given in the text by putting a’=s?/3 
and b’?=30?/2. Both the formulae for p, are valid in any 


coordinate system but for the special one of the principal 
axes the invariant b’? becomes 


b?= $[(arz— ayy)? + (ty y— zz)? + (azz—@rz)* ]. (4a) 
b’ is called the optical anisotropy. 
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CHLORINE VIBRATIONS 





IN DICHLOROBENZENES 


TABLE I. Transformation coefficients for the Raman tensor azy. 














Orr ayy Qzz ary Qyz ” axe zr Qyz zy 
att cos? » sin? @ 0 cos gsin ¢ sin g COs ¢ 0 0 0 0 
Onn sin? 9 cos? 9 0 —sin gcos ¢ —cos gsin ¢ 0 0 0 0 
age 1 0 0 0 0 0 
Qty —cos g sin ¢ sin g cos ¢ 0 cos? ¢ —sin? 9 0 0 0 0 
Ont —sin gcos ¢ cos g sin ¢ 0 —sin? ¢ cos? 9 0 0 0 0 
ary 0 0 0 0 0 cos ¢ 0 sin ¢ 0 
are 0 0 0 0 0 0 cos ¢ 0 sin ¢ 
Ont 0 0 0 0 0 —sin ¢ 0 cos ¢ 0 
Qty 0 0 0 0 0 —sin 9 0 cos ¢ 








two deformation vibrations in the molecular 
plane and the two deformation vibrations out of 
the plane. The results will be compared with 
experimental data. 

The vibrations with which we are concerned 
here are represented schematically in Fig. 1 to- 
gether with their symmetry symbols.* In each 
pair of vibrations the two chlorine atoms vibrate 
in phase for one vibration and out of phase for 
the other one (called A and B vibrations respec- 
tively in Fig. 1). If we assume that in first ap- 
proximation the Raman tensor for these vibra- 
tions is determined by the change of polariza- 
bility in the C—Cl bonds alone and moreover 
that the contributions from the two bonds are 
additive, then we can construct the Raman ten- 
sors by either adding or subtracting the single 
C—Cl contributions with reference to a chosen 
system of axes rigidly connected with the mole- 
cule. This system with axes &, 7, ¢ is indicated in 
Fig. 1. € and 7 lie in the molecular plane, ¢ is 
perpendicular to it. If the Raman tensor of a 
particular C—Cl bond has the six components 
az, With the x axis in the C—Cl direction, the 
y axis perpendicular to it in the plane and the 
z axis orthogonal to the plane then we may ex- 
press this tensor in the coordinate system &, n, ¢ 
by making the usual tensor transformation 


Oey = x Qzy COS (x£) cos (yn). (5) 


From the special choice of the six axes it follows 
that all terms in the sum containing factors 
cos (x), cos (¢y), cos (£z), and cos (nz) are zero. 


3 The symmetry symbols used are listed for instance in 
H. Sponer and E. Teller, Rev. Mod. Phys. 13, 76 (1941). 
See also A. Langseth and R. C. Lord, Kgl. Danske Vid. 
Sels. Math.-fys. Medd. 16, No. 6 (1938) for the classifica- 
tion of the vibrations of the deuterobenzenes according to 
symmetry and type, and H. Sponer and J. S. Kirby-Smith, 
reference 1, for the vibrations of the chlorobenzenes. 


If yg is the angle (xt) we may express all the 
coefficients occurring in the transformations (5) 
as functions of ¢ as indicated in Table I. 


I : . 
The Raman tensor a,, belonging to the dis- 
placement of the chlorine atom Cl; becomes, with 


the help of the angle of transformation g, ees 
in the reference system £, 7. Similarly the angle 
of transformation is —¢ for the chlorine atom 
Clz;. To indicate this, the corresponding tensors 


I II ‘ id 
are called a:,(¢) and a:,(—¢). Their addition and 
subtraction gives the resulting Raman tensors 
for each pair of vibrations A, B. 


I II 
A= ae (¢) + aey( —¢), 
I IT (6) 
B=ay,(¢) Faiq(—¢). 
The + and — signs in (6) have to be determined 
for each of the three pairs separately (see below). 
The intensity ratio of the two vibrations v4 
and vg can be calculated from Eq. (3) and the 
depolarization factors ps and pz from Eq. (4), 
after A? and B? have been determined. 


A. The Chlorine Valence Vibrations 


We assume, as is customary, that for small 
amplitudes the change of polarizability varies 
linearly with nuclear displacements. In our case 
this change for one C—Cl group may be taken 
as the one which occurs during the totally sym- 
metrical chlorine vibration in CsH;Cl. If we make 
the simplifying assumption that the Raman ten- 
sor az, corresponding to this vibration has cylin- 
drical symmetry about the C—Cl axis (which we 
take as x axis), then the tensor has the form 


ae +8 8 aan © Of 
Aay=|0 ay O|/=|0 ay O}. (7) 
0 O azz 9 O80 @& 








We can determine the ratio a;:/a,, from the 
experimentally known‘ depolarization factor 
pn=0.13. From (4) and (4a) follows then 





with 
$?=3 (22+ 2ayy)? and }?= §(Qz2— yy)”. 


Since we are interested in the ratio s?/b? only, 
i.e., in relative intensity values, we may arbi- 
trarily fix one constant, for example s?=3. With 
this value we obtain for the intensity using 
Eq. (3) 

Tmono™~ 10s?+ 136? = 30+10=40 


and, if we take the positive root of s?/b?, we have 
Ge2=1.716; ayy=asg=0.642. 


This choice does not mean a restriction of 
generality in the final results given below because 
the intensity and depolarization factors always 
depend on s?/b? only, and the negative root leads 
to the same result. The shape of the polarizability 
ellipsoid, however, depends on the sign. The 
positive sign, yielding a,,>ay,=a-.- is physically 
more plausible since it corresponds to an ellipsoid 
with the major axis in the direction of the C—Cl 
bond. 

In order to find the Raman tensors belonging 
to the two valence vibrations A and B of the 


dichlorobenzenes (Fig. 1a), we calculate ies and 


Psi in the coordinate system &, 7 from (7) with the 
help of Table I. We can see from Fig. 1a that the 
Raman tensor A belonging to the totally sym- 
metrical valence vibration, in which the chlorine 
atoms are in phase, will be given by adding 


T IT 
a:,(y) and a:,(—) and the tensor B (chlorine 
atoms out of phase) by subtracting them. 
* (6) becomes here 


I II 
A =az,(~)+az,(—¢), 
I II 
B=ay(¢) —ae(—¢), 
with 
Az:=2(az2 cos* g+a,, sin? ¢), 


A m= 2(azz Sin? g+ay, Cos? ¢), 
Ax =2az:, 


‘L. Simons, Soc. Sci. Fennica, Comm. Phys.-Math. 6, 
No. 13 (1932). 
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An=Ay=Ag=Ay=Ay=A;,=0, 
By =By=By=0, 

By, = By = —2(ar2— ayy) Sin g COS ¢, 
By = By: = By = B;,=0. 


With the values of azz, ayy, a2 determined 
above we find for all three dichlorobenzenes s 4? 
= (2s)?=12; ss*=0. 

The squares of the tensors are given by 


A?=A4[ 3.36(cost g+sin‘ ¢) 
+4.4 sin? yg cos? 9+0.41], (8) 
B?=4X2.32 sin? ¢ cos? ¢. 
Thus 
b4?=A?—s,4?=4[ 3.36(cost o+sin! ¢) 
+4.4 sin? ¢ cos? g—2.59 ], 
bp? = B*=4X 2.32 sin? ¢ cos? ¢. 


The intensities J, and J and their ratios can be 
found with the help of (3) and the depolarization 
factors pa from (4). Since sz?=0, it follows 
pp=6/7 which means that the non-totally sym- 
metrical vibrations are completely depolarized as 
expected. 

The intensity ratios and depolarization factors 
for the vibrations of the three individual di- 
chlorobenzenes can be calculated with the above 
expressions by putting g=90°, 60°, 30° respec- 
tively for the para- meta-, and ortho-compound, 
(comp. Fig. 1(a)). The results are summarized in 
Table II together with the data on the totally 
symmetrical valence vibration of monochloro- 
benzene. The discussion of this table will be given 
in a later section. 


B. The Chlorine Deformation Vibrations 


We assume that the main effect of the chlorine 
deformation vibrations consists in a rotation of 
the polarizability tensors a connected with the 
C—Cl bonds. We neglect the distortion of a due 
to changes in the C—Cl distance. This is analo- 
gous to the procedure for the valence vibrations 
where we considered the change in distance only 
and neglected the angular changes. However, it 
might lead to less reliable results because the 
distortion of the carbon-ring might give a rela- 
tively larger contribution to the Raman tensors 
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TABLE IT. Calculated relative intensities and depolariza- 
tion factors for the chlorine valence vibrations of the 
dichlorobenzenes. 








Chloro- 





benzene IA IB Ta/IB PA PB 
p-—di 4x40 0 a 0.13 —- 
¢g=90° 

m—di 

y= 60° 4 34.4 4X 5.56 6.1 0.063 6/7 
o—di 

g=30° 4X 34.4 4X 5.65 6.1 0.063 6/7 
mono 40 —_ — 6:13" — 








* Taken from experimental data (reference 4) to adjust the constants. 


of the deformation vibrations than the angular 
changes to the valence vibrations. 

A deformation vibration consists in a rotation 
of the C—Cl bond (axes X, Y, Z) against a 
system of axes x, y, z in which the C—Cl bond is 
at rest. We have to calculate first the Raman 


tensor ie (change of polarizability) produced by 
such a rotation of a single C—Cl group. The effect 
of the two C—Cl groups in the dichlorobenzenes 
can then be composed with the help of (6) in the 
same manner as for the valence vibrations by 
transformation into a coordinate system &, 7, ¢ 
fixed in the molecule. 

For the deformation vibrations in the plane the 
rotation of one C—Cl group will take place in 
the x, y plane (8 =angle of rotation), the z and Z 
axes being coincident. The coordinate systems 
used for the calculation of the Raman tensors be- 
longing to these vibrations are shown in Fig. 1(b). 
The calculation for the non-planar deformation 
vibrations is analogous; there the rotation of one 
C—Cl group takes place in the x, z plane per- 
pendicular to the plane of the molecule and the y 
and Y axes coincide. (See Fig. 1(c).) 

The polarizability tensor axy rigidly connected 
with one oscillating C—Cl bond (diagonal terms 
only) can be expressed in the system x, y, 2 
(C—Cl bond at rest) by again using Table I. The 
coefficients will be functions of the angle of rota- 
tion #. Since the rotation of the C—Cl axis out 
of the position of equilibrium will be small, the 
angle & will be small and may be written as 
3=e€cos 2mvit where v; is the frequency of the 
vibration. We can then develop the tensor az, 
in e and only keep linear terms in ¢ (sin 89, 


cos 31). Corresponding to the second row in Eq. 
(1) the terms of zeroth order give the Rayleigh 
scattering and the terms of first order the com- 
ponents of the Raman tensors. 

For a C—Cl bond rotating in the x, y plane we 
have 


T - T 
Arr = AXX cos” d+ayy sin? d~axx, 


—_— = T e 
Qyy, = axx sin? d+ayy cos? d~ayy, 
T 
Azz=AzZZ; Or: = Azz = Ayz = A2,=0, 


T T = 
Ary = Ayr = (—axx+ayy) sin ? cos 3 
T T 
= (—ayx+tayy)e cos 2mp;jt. 


For the angle of rotation —#, az,= ay: changes its 
sign; the other components remain the same. 

In the system where the molecule is at rest, 
the Raman tensor belonging to one C—Cl group 
vibrating in the plane will have the components 


R R R R R R R 
Arr = Ayy = Azz = Agz = Azz = Ayz = Ary = 0, 


(9) 


R T 
Ory = Ayz=(—axxtayy)e=C. 

A similar calculation for the rotation of a 
C—Cl bond in the x, z plane leads to the follow- 
ing components of the corresponding Raman 
tensor 

R R R R R R R 
Ope = Ayy = 22 = Ary = Ayz = Ayz = Az, =0, 


R R T a ’ ’ 
Orz=Q22= (axx—azz)e =c’. 


(10) 


cand c’ in (9) and (10) are certain unknown con- 
stants in which the intensities of the vibrations 
can be expressed. 

The Raman tensors of the deformation vibra- 
tions of the dichlorobenzenes can now be easily 


obtained. a and een denote again the Raman 
tensors belonging to the two C—Cl bonds in the 
coordinate system £, 7, ¢. These tensors can be 
found by transformation with the angle ¢ of 
tensors of the form (9) (for the planar deforma- 
tion vibrations) or (10) (for the non-planar vibra- 
tions) using again Table I. The dependence on the 
angle can be read from Figs. 1(b) and 1(c). We see 
that Eqs. (6) become for the planar vibrations 


T II 
A =az,(v, ¢) tae(—9, —¢) 
I Tl 
=az,(0, ¢) —az, (8, —¢), (lla) 
I IT 
B=az,(d, ¢) +az,(8, —¢), 
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TABLE III. Calculated intensities and depolarization 
factors for the planar chlorine deformation vibrations of the 
dichlorobenzenes. 








Chloro- 





benzene IA IB TA/IB pA PB 
p-—di 
¢=90° 0 13X 8c? 0 — 6/7 
m—di 
¢=60° 13 X 6c? 13X2¢c 3 6/7 6/7 
o—di 
g=30° 13 X 6c? 13X 2c? 3 6/7 6/7 
mono — 13X 2c? - — 6/7 








and for the non-planar vibrations 
I a 
A = en(B, ¢) +ae,(d, Nad ¢), 
I II 
B= aen(B, ¢) +az(—9, —¢) 
I IT 
=az(3, y) —ae,(3, —¢). 


For the pair of vibrations in the plane the com- 
ponents of the Raman tensors become 


(11b) 


Au=a,{2 cos g sin g—2 cas (—¢) sin (—¢) ] 
=c[4 cos ¢gsin ¢ ], 

hata, —-2 cos ¢ sin g+2cos (—¢) sin(—¢) ] 
=c[{—4cos gsin ¢ ], 

Ay, =Ay: 
=«a,,[ cos? g—sin? g—cos* (— g)+sin® (—¢) ] 
=0, 

Ay =Ag=Ay=Ay=A;,=0, 

Bye= By = By = By = Bye = By = By, = 0, 

Bey = By = 0242[ cos? g—sin? ¢ | 
=c2[cos? g—sin? g]. 

As before, we obtain a measure for the inten- 
sity with the help of (3) after determining the 
squares of these tensors 

A*=c? X32 cos? ¢ sin? ¢, 
B?=c* X8(cos* g—sin? ¢)?. 

Since s42=5,7=0, wehavel,=13A?, Jp =13B?, 
and [4/Ip=A?/B?*. From Eq. (4) we have pa=pz 
=6/7 which means that all the deformation vi- 
brations in the plane are completely depolarized. 


For the deformation vibrations out of the plane 
the components of the Raman tensors become 


(12) 








AND H. SPONER 


with the help of (10), (11b), and Table I 
Az=Am=Ayz=Apy=Aye=0, 
Ag = Aye=az-[c0s gt+cos (—¢) ]=c'(2 cos ¢), 


Ag=Ap,= ar, —sin g—sin (— ¢) ]=0, 
By = By = By = Bey = By = By = B;, =0, 


By = Bry =ae.[ —sin gt+sin (— yg) ]=c’(—2sin ¢). 


Since s4?=s,*=0, the intensities are again di- 
rectly proportional to the squares 


I,4=13A?=13c’? X8 cos? ¢, 


. (13) 
Ip=13B?=13c'?X8 sin? ¢. 


The intensity ratio is given by cos? g/sin? g and 
pa=pa=6/7. 

The intensity ratios for the pairs of deforma- 
tion vibrations of the three dichlorobenzenes can 
be obtained by putting ¢=90°, 60°, 30° respec- 


TaBLE IV. Calculated intensities and depolarization 
factors for the non-planar chlorine deformation vibrations 
of the dichlorobenzenes. 








Chloro- 
benzene Ia IB Ia/IB pA PB 





p—di 


¢=90° 0 13 X 8c” 0 _ 6/7 


m—di 


¢=60° 13 X 2c” 13X 6c” 1/3 6/7 6/7 


o—di 


¢= 30° 13 X 6c” 13X 2c” 3 6/7 6/7 


mono 13X2c’2 — — 6/7 — 








tively. The results are summarized in Tables III 
and IV. Their discussion and the comparison 
with the experimental data will be given in the 
next section. 


DISCUSSION AND COMPARISON WITH 
EXPERIMENTAL DATA 


To facilitate the discussion of the results of the 
previous section in relation to existing experi- 
mental data, all observed low frequency Raman 
lines (up to the strong 1000 cm™ region which 
corresponds to the totally symmetrical ring vi- 
brations) of the dichlorobenzenes have been col- 
lected in Table V, together with the estimated in- 
tensities and observed depolarization factors. All 
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TABLE V. Raman spectra of dichlorobenzenes. 








, Raman lines 
Swaine and 


Depolarization factors 
Sponer and 











Murray Paulsen Kirby-Smith Paulsen 
Substance (Ref. 5) (Ref. 2) (Ref, 1) (Ref. 2) Cale. Assignment 
219(00) 
p—C.H.Cl, 232(00) 
302(4) 298(65) 0.8 0.73 0.86 Bog (C—Cl def.) 
¢=90° 333(8) 330(134) 0.5 0.21 Qig (€g* 
355(1) 0.86 Bay (C—Cl def.) 
479(00?) 
509(00?) 
627(5) 631(36) 0.8 0.78 Bog (€g*) 
710(00) 
748(105) 748(104) 0.2 0.06 0.13 aig (C—Cl valence) 
794(00) , 
812(1) 
Dadieu, 
Pongratz, and 
Kohlrausch 
(Ref. 7) 
178(3) 176(5) 0.9 0.86 a2 a1 
m—C.H,Cle 202(3) 199(8)d) 0.9 0.86 Be or ay (3 Or a2 
216(2) 210(8d) 0.86 a, Or Bo \ ae or Bi 
¢=60° 366(1) 360(3) 
399(4d) 397(8) 0.10 a1 (€,*) 
428(2) 428(4) 0.6 Bi (€*) 
530(0) 
666(4) 662(6) 0.25 0.06 a: (C—Cl valence) 
154(10) 152(10b) 0.8 0.86 Bo a1 
o—C,H,Cl. 203(3) 199(4) 0.8 0.86 a2 Bi or Be 
239(1b) 240(5d) 0.8 0.86 a Be or Bi 
g=30° 330(0b) 336(3) 
430(1) 427(5) 0.8 Bi (€,*) 
469(1) 469(5) 
483(2) 486(6) 0.45 a1 (€,*) 
658(5) 660(8) 0.35 0.06 a, (C—Cl valence) 
756(0b) 747(1/2) 
860(1d) 








measurements refer to the liquid. Swaine and 
Murray’s® measurements of the lines and Sponer 
and Kirby-Smith’s! measurements of the depolar- 
ization factors are given for all three substances. 
For para-dichlorobenzene a careful study with 
improved technique has been published by Paul- 
sen.? We have reproduced here all his lines, also 
the very weak and doubtful ones, together with 
his depolarization factors which at present repre- 
sent the most reliable values. The calculated de- 
polarization factors have been taken from Tables 
II-IV of the preceding section. The assign- 
ments are partly the same as in reference 1 and 
are partly suggested from our theoretical results. 

The discussion will be carried through accord- 
ing to types of vibrations as in the previous 
section. 


5 J. W. Swaine and J. W. Murray, J. Chem. Phys. 1, 512 
(1933). 


A. The Chlorine Valence Vibrations 


As mentioned in the Introduction (see also 
Fig. 1(a)), there are two types of valence vibra- 
tions in each dichlorobenzene which are dis- 
tinguished by motions of the two chlorine atoms 
in phase and out of phase (generally denoted by 
A and B in this paper). As pointed out in refer- 
ence 1, they originate from a linear combination 
of benzene hydrogen vibrations in the 3000 cm~ 
region, when replacing two hydrogen by two 
chlorine atoms, and they belong to the following 
symmetry classes: ai, (para) from aig, €,* in 
benzene, 81. (para) from Bix, €.~; a, (meta) from 
O19, Biu, €9*, €u-; Bi (meta) from e,t, €.~; and a; 
(ortho) from ayy, €g*, €u~; 81 (ortho) from Byu, €,*, 
¢.-. The totally symmetrical vibration has been 
observed with certainty in all three cases in the 
lines 748 cm (a1, para), 666 cm (a; meta), and 
658 cm (a; ortho), as has been noticed by many 
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authors. The intensity estimates of these lines 
given in Table V agree with our calculations, as 
far as the ratios between the lines of the various 
compounds are concerned. The lines in the meta- 
and ortho-compounds are of about the same in- 
tensity while the corresponding line in para- 
dichlorobenzene is stronger. All three lines are 
highly polarized. The depolarization factor of a1, 
(para) should theoretically coincide with p=0.13 
of a; in CsH;Cl which we have used as fixed 
value.* Sponer and Kirby-Smith’s value is slightly 
higher and Paulsen’s more accurate figure is 
lower than the one calculated by our rough 
method so that the agreement can be called satis- 
factory. The p’s for the a;-line in the meta- and 
ortho-compound should theoretically have the 
same value, namely about $ of that for a1, (para). 
This result is not confirmed experimentally. Al- 
though, according to the experimental method, 
the measured p-values represent upper limits and 
more accurate values will probably be noticeably 
lower (analogous to Paulsen’s p=0.06 of ar, 
(para) instead of 0.2), it is improbable that 
they will be below the p-value for the a1, para- 
frequency. The assumption of cylindrical sym- 
metry for the change of polarizability tensor for 
one vibrating C—Cl group is conceivably not a 
good approximation for this vibration. However, 
the assumption cannot be very wrong and it 
seems to us that a slight deviation from it does 
not alter the results in the desired direction so 
that other effects like non-additivity might be 
responsible. 

An assignment for the unsymmetrical valence 
vibration which must have a similar frequency 
value has never been reported with certainty. For 
the para-compound the vibration §;, is forbidden 
in the Raman effect but may occur in the infra- 
red. Lecomte® assigns the symmetry #1, to the 
weak infra-red line 620 cm~!. This interpretation 
is somewhat doubtful. The fact that this line 
changes its frequency value somewhat with the 
substituents is in favor of his assignment; its low 
value is against it since the unsymmetrical val- 
ence vibration is frequently larger than the sym- 


* The intensity of the Cl-vibration a; in CgH;Cl of fre- 
quency value 702 cm™ is slightly weaker than aig in 
p-CsH,Cls, although the intensity ratio should theoretically 
be about }. However, a comparison is difficult since the 
intensity estimates were made by different authors. 

6 J. Lecomte, J. de phys. et rad. 9, 13 (1938). 
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metrical one. Besides the infra-red band just 
mentioned, Lecomte observed quite a number of 
bands in the region in which we are interested, 
namely a strong band at 820 cm and weak ones 
at 650, 680, 720, 750, and 780 cm~. He inter- 
preted the bands at 820 and 720 cm~ as coming 
from the non-planar hydrogen vibrations e,*+ 
(about 850 cm~') and ag, (671 cm infra-red) 
in benzene. Both of them lead to 82, vibrations 
in the para-compound which are allowed in the 
infra-red. However, this assignment cannot be 
correct, as one of these two 82, vibrations be- 
comes a non-planar chlorine deformation vibra- 
tion which must be considerably lower in fre- 
quency. We would prefer to assign the weak 620 
cm! in the infra-red to the same vibration as the 
medium strong Raman line at 631 cm™, namely 
to the B2, (coming from the e,+ in benzene), which 
is forbidden in the infra-red but might appear 
weakly in the liquid on account of intermolecular 
forces. The strong infra-red line 820 cm™ could 
be assigned to the valence vibration 8,; this vi- 
bration might be accompanied by a sufficiently 
strong variation of the electric moment so that 
the corresponding infra-red band could be fairly 
intense. Paulsen’s very faint line at 812 cm 
might be caused by the same vibration which is 
forbidden in the Raman effect but might appear 
weakly because of intermolecular forces. The 
band 720 cm“ (infra-red) and the extremely faint 
Raman line at 710 cm™ reported by Swaine and 
Murray might be due to the other 62, vibration 
as suggested by Lecomte. For his tentative as- 
signment of the weak bands 680 and 750 cm 
to non-planar 83, vibrations (82, and e,~ H-vibra- 
tions in benzene), forbidden in the infra-red, the 
same is true as for the 62, vibrations; one of them 
should become a Cl-bending vibration with 
smaller frequency. On the whole, it seems prema- 
ture to offer definite assignments for all weak 
lines, particularly because selection rules may 
break down in the liquid and thus may cause 
forbidden lines to appear in the Raman effect or 
infra-red respectively. 

In the meta- and ortho-compounds the unsym- 
metrical C—Cl valence vibration has the sym- 
metry §; and it is allowed although the corre- 
sponding line is much weaker than the a. It is 
completely depolarized in the Raman effect. It 
should have a frequency value similar to the line 
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of the totally symmetrical valence vibration. We 
see from Table V that no lines of reasonable in- 
tensity are available for this interpretation. We 
already know from Table II that the intensity 
ratio between the two valence vibrations should 
be about 6, i.e., the symmetrical valence vibra- 
tion should occur about 6 times stronger than the 
unsymmetrical one. This places the expected 
Raman lines at once among the very weak ob- 
served lines. The only lines which occur in the 
right region are the weak lines at 756 and 860 
cm~! in ortho- and possibly 530 cm in meta- 
dichlorobenzene. However, a number of weak and 
stronger bands have been reported in the infra- 
red. Lecomte gives tentative explanations for all 
of them. According to him the lines at 756 and 
860 cm, both of which appear also in the infra- 
red, result from the H-vibrations e,~ and e,*+ in 
benzene, respectively. The interpretation is pos- 
sibly correct. Lecomte places the 6; chlorine va- 
lence vibration at about 820 cm in both the 
meta- and ortho-compounds. Since this fits in 
well with our interpretation of the same line in 
the para-compound, we consider this as probable. 
From our computations we can only conclude 
that the unsymmetrical halogen valence vibra- 
tions will not show up at all or as very weak 
Raman lines although they are allowed by selec- 
tion rules. 


B. The Chlorine Deformation Vibrations 


The two groups of deformation vibrations of 
the three dichlorobenzenes (in the molecular 
plane and out of it) belong to the following sym- 
metry classes: Planar vibrations have the sym- 
metries 83, and 82, for the para- and a and #; for 
the meta- and ortho-compounds. Non-planar vi- 
brations have the symmetries 62, and 83, (para) 
and $2 and a2 (meta and ortho). The two types 
are denoted with A and B in this paper and are 
shown in Fig. 1(b) and 1(c). The monochloroben- 
zene has only one type of chlorine vibration in 
each group of symmetry, @; (planar) and B:2 (non- 
planar). It is unessential for our discussion 
from which benzene vibrations these vibrations 
originate. 

As can be seen from Tables III and IV we ex- 
pect the following results from our computations 
regarding the appearance of these vibrations: The 
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vibrations 83, and B2, of p-dichlorobenzene are 
forbidden in the Raman effect. For the planar 
vibrations, the a;’s in m- and o-dichlorobenzene 
should appear with the same intensity and should 
be three times as strong as the 6; vibrations. The 
82, para-vibration should be four times as strong 
as these §;-vibrations. For the non-planar vibra- 
tions, we expect the az (meta) to be of the same 
intensity as the ®2 (ortho) and three times as 
strong as the 82 (meta) and ae (ortho). The #3, 
para-vibration is, according to calculations, 
slightly stronger than a2 (meta). 

In our approximation all the deformation vi- 
brations give completely depolarized Raman 
lines, including the totally symmetrical ones, so 
that we cannot distinguish between them in this 
respect. The reason for this is, of course, our ne- 
glecting all deformations of the change of polari- 
zability tensors and considering their rotation 
only. Such a computation leads to a zero trace of 
the tensor. A better approximation which takes 
the deformations into account would lead to a 
small trace value for the totally symmetrical 
vibrations so that their depolarization factors 
would be slightly less than 6/7. 

We turn now to a comparison of these theo- 
retical results with the observed low frequency 
Raman lines of the dichlorobenzenes. Para- 
dichlorobenzene has a strong and highly polarized 
line at 330 cm which has been interpreted as 
caused by a symmetrical carbon vibration (re- 
sulting from ¢,+ in benzene). A line of medium 
intensity was found at 298 cm by all observers. 
It occurs also in the ultraviolet spectrum and 
might well correspond to the $2, vibration. A 
weaker line has been reported at 355 cm™! by 
Swaine and Murray but not by Paulsen. It repre- 
sents perhaps the 83, vibration. The extremely 
faint lines which Paulsen reports at 219,232,479, 
and 509 cm must correspond to vibrations for- 
bidden in the Raman effect. Besides the two Cl- 
deformation vibrations of symmetry 83. and Bou, 
carbon frequencies of symmetry Bo. and ay (€u* 
in benzene) might appear weakly. Furthermore, 
besides the already mentioned carbon vibration 
B2,(€,*+) of probable frequency 631 cm™, another 
allowed carbon vibration of symmetry B39 (G2, in 
benzene) should occur in this region. The infra- 
red band at 540 cm of medium intensity re- 
ported by Lecomte might correspond to Paulsen’s 
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509 cm and could be any 8, vibration since they 
are allowed in the infra-red. 

The correlation is more difficult for the m- and 
o-dichlorobenzene since all vibrations are allowed 
in the Raman effect. Also, the intensity estimates 
are less accurate than for the p-compound. The 
essential question is, which type of bending vibra- 
tion belongs to the lowest frequency value. The 
values lie close enough to make either choice pos- 
sible. We have pointed out before that the calcu- 
lated intensity ratios are the only indication 
which the theory gives us: a; : 8: =3 : 1 (mando) 
and a2(m) : B2(m)=a2(m) : a2(o0)=3 : 1, where 
a2(m) has the same intensity as 62(0). Consulting 
Table V we may first rule out the frequencies 399 
cm in m—CgHyCle and 483 cm in o—CgHuCle 
which belong to the symmetrical carbon vibra- 
tion a; resulting from the 606 e,+ vibration in 
benzene. The non-symmetrical component lies in 
the same frequency region and is probably the 
428 in meta and the 469 or 430 cm in ortho. 
This leaves us with the lines 178(3), 202(3), 
216(2), 366(1), 530(0) for m—C,H,Cl. and 
154(10), 203(3), 239(1b), 330(0b), and 430(1) or 
469(1) in o—C,.H,4Cle, from which we have to 
select the lines belonging to the deformation vi- 
brations. We see immediately that it is impos- 
sible to fulfill the theoretical intensity relations. 
Several different sets give more or less equally 
good fits depending on whether one uses the in- 
tensity estimates by Swaine and Murray® or by 
Dadieu, Pongratz, and Kohlrausch’ throughout 
for consistent comparisons. Each set has its own 
advantages and disadvantages. As long as we do 
not have more precise intensity measurements, 
it seems useless to prefer one of these correlation 
sets. We must try to take recourse to other 
sources. The ultraviolet absorption spectra of all 
three dichlorobenzenes® have been studied. From 
the occurrence of the low frequencies in these 
spectra one can draw certain conclusions about 
the symmetries of the vibrations concerned. Com- 

7A. Dadieu, A. Pongratz, and K. W. F. Kohlrausch, 


Monatsh. f. Chemie 61, 426 (1932). 
8H. Sponer, Rev. Mod. Phys. 14, 224 (1942). 
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bining all the information available we would 
believe that the sets, in which we identify the line 
154 cm™ (0) with 6B. and 178 (m) with a2 are 
more probable than sets with 154 (0) and 178 (m) 
as a;. Probable assignments, giving more weight 
to the more recent measurements, are indicated 
in the last column of Table V. We refrain from 
more definite assignments until accurate intensity 
measurements are available. We hope that then 
such simple theoretical calculations as presented 
here may be of help for the interpretation of 
Raman spectra of these and similar molecules. 

After conclusion of our calculations, a paper 
by Wolkenstein® came to our attention which has 
some general bearing on the question discussed 
here. Wolkenstein uses the same principle, but in 
a more elaborate way, of calculating Raman in- 
tensity and depolarization factors for simple 
molecules. He also assumes additivity of the 
Raman tensors of individual bonds, but he uses 
a complete system of force constants to deter- 
mine the actual deformations and displacements 
in the various normal vibrations. This is, of 
course, a more correct method, if it can be carried 
through. It seems, however, that the results ob- 
tained are not superior to those of our simplified 
procedure. Wolkenstein treats the chlorine and 
bromine vibrations of CH3;Cl, CH3Br, CHeCle, 
and CH:2Brz. As an illustration for the applica- 
bility of our method, we calculated the depolari- 
zation factors and relative intensities of the two 
C—Cl valence vibrations of CH2Cl: using the 
depolarization factor of the C—Cl valence vi- 
brations of CH;Cl, and obtained results which 
are in as good agreement with observations as 
Wolkenstein’s. 

In conclusion we should like to express our 
thanks to Dr. E. Teller who stimulated the theo- 
retical treatment of the problem presented here 
and to Dr.-L. W. Nordheim for helpful discus- 
sions. We also wish to acknowledge a grant-in-aid 
from the Duke University Research Fund for 
carrying out this research. 


( 9M. W. Wolkenstein, J. Phys. Acad. Sci. U.S.S.R. 5, 185 
1941). 








i- 








THE JOURNAL OF CHEMICAL PHYSICS 





VOLUME 11, NUMBER 6 JUNE, 1943 


Optical Sensitizing of Silver Halides by Dyes 
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Absorption and Constitution* 
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A comparison of increasingly halogenated fluoresceins 
shows that, with increasing weight of halogen introduced, 
the principal visible absorption moves steadily to longer 
waves, as also the first ultraviolet band. At the same time 
the acid dissociation constants of the phenolic and car- 
boxylic groups tend to become indistinguishable, while the 
solubilities of the silver compounds decline. The various 
structures involved in the resonance of the alkaline, 
hydrogen, and silver compounds are discussed. Measure- 
ments of the solubility of silver erythrosinate show this to 


vary with ‘aging’ and particle size. The solubility of 
silver erythrosinate on adsorption to silver bromide is 
apparently lower than that of the pure compound. The 
adsorption is oriented, an edge-on or end-on monolayer 
being formed. Although not irreversible, there is a peculiar 
deviation from the Langmuir isotherm. The optical 
sensitizing appears to occur by a mechanism similar to 
that with basic dyes. Initially, transfer of an exciton to 
the silver halide seems indicated, rather than of an electron. 





1. THE ABSORPTION SPECTRA OF FLUORESCEIN 
DYES IN RELATION TO CONSTITUTION 
AND SOLUTION 


ERTAIN fluorescein dyes, viz., eosin and 
erythrosin, are the oldest of the practical 
color sensitizers for silver halide emulsions.! 
Although they have been displaced almost 
entirely by cyanine (polymethine) dyes, even in 
their own region of green-yellow (orthochromatic) 
sensitizing, they exhibit certain characteristics 
of definite importance for the basic theory of 
optical sensitizing. Some essentials of the mo- 
lecular constitution of dyes of this category have 
been excellently described by G. Schwarz.? 
The archetype of the substances is xanthene 


O 


H: 
(acid form) 
in which, however, the O-atom can be replaced 
by such bivalents as S, Se, CH. ---. Of more 
immediate importance is the genotype fluorescein 
which is (water) soluble by ionization with 
alkali of the —OH group in the xanthene 
skeleton, and of the —COOH group in the 
pendent phenyl nucleus. Also the xanthene 


*Communication No. 914 from the Kodak Research 
Laboratories. 

1Cf. S. E. Sheppard, Atti. X. Congr. Internaz. Chim. 
Rome 1, 234 (1938). 

2G. Schwarz, Sci. Ind. Phot. 8, (2) 97 (1937). 
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skeleton has become an extended resonance sys- 
tem because of the =O : —O~ groups. For all 


O O O- 


of this, fluorescein is quite incompetent as an 
optical sensitizer for the silver halides. This 
property appears, with strengthening of the 
resonance system, by the introduction of halogen 
atoms in the xanthene yg —_— 


ua O A/ 
O 
the tetrabrom- _. 
fluorescein COOH 
i 
Eosin 
O ; & 2 OH 
\ ff 
and the I I 
tetraiodo- ae 
fluorescein COOH 
Yd 
Erythrosin 
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NATE 
Scm TUGE 


o pH*=4.0 
x pHe3.0 
4& pH *12.0 


DENSITY 





500 
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Fic. 1. 


are two important sensitizers of this group. It is 
interesting to observe what occurs by intro- 
duction of such electronegative atoms in the 
xanthene frame. Table I shows that the spectral 
absorption is only slightly shifted toward the 
longer waves, and it is evident that the funda- 
mental resonance system is little affected, be- 
cause the ratio of the wave-lengths (or fre- 
quencies) of the visible to the first ultraviolet 
band is substantially unchanged. This first 
ultraviolet band probably represents the second 
excitation level. 


Acid Strength 


The fluoresceins are acid dyes, and have two 
potentially acidic groups, viz., the phenolic —-OH 
of the xanthene frame, and the —COOH group 
in the dependent phenyl nucleus. Attempts were 
made to determine the corresponding dissociation 
constants (see Table II) by electrometric titra- 
tion of the alkaline solutions with acid. Only in 
the case of fluorescein itself were definite values 
obtained for the two groups.* In the other cases 
values were obtained of one pK value; this 


* Measurements by R. E. Stauffer of these Laboratories. 
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would indicate that in these the pK values were 
too close to be separated. In fact, the solubilities 
of the later acids in water are so low that the 
pIK values given may be dictated by the solu- 
bilities, and not express the actual dissociation 
values. The acidification is accompanied by a 
change in the spectrum, which, in water, is 
chiefly that occasioned by withdrawal of the 
fluoresceinate ion from solution. While, for 
example, the erythrosin acid has (in alcohol) an 
absorption differing considerably from that of 
the ion, it can contribute little to the spectrum 
in water because of its low solubility. Figure 1 
shows the nature of the change as the fH is 
lowered. If the pH is maintained sufficiently 
high (cf. Fig. 2), erythrosin in alkaline solution 
—as sodium erythrosinate—follows Beer’s law 
up to fairly high concentrations, when secondary 
“salt” effects show themselves. These have been 
studied in the case of eosin (tetrabromfluo- 
rescein) by G. Kortum* who attributed devia- 
tions from Beer’s law to the formation of ion 
aggregates. However, eosin and erythrosin show 
no such change in spectral distribution (intensi- 
fication of $-band relative to a-band) with 
increased concentration as is found with thiazine 
and cyanine dyes, and which is attributed to 
dimerization.‘ One of us has indicated elsewhere‘ 
that it is difficult for the pendant phenyl! nucleus 
to be coplanar with the xanthene frame, and 
that this would keep the erythrosinate ions from 
parallel apposition at less than 4—5A interval 
between them. This distance would probably 
preclude the spectral anomalies shown in dimer- 
ization. The ‘‘fluoresceinate” ion can have the 
two structures shown below: 


O O om O- O oO 
COO COW 

\ \ O S44 O 
A | fo 


a” i 
O- O- 


These would have nearly equal energies, and 
constitute a resonance system accounting for the 
absorption band in the visible region. The 


3G. Kortum, Zeits. f. physik. Chemie 34, 255. (1936). 
*Cf. S. E. Sheppard, Rev. Mod. Phys. 14, 305 (1942). 
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introduction of electronegative groups in the 
xanthene nucleus would produce secondary elec- 
tromeric effects’ somewhat strengthening the 
resonance but not essentially changing the sys- 
tem. With the addition of protons to form the 
(undissociated) acids certain complications oc- 
cur. In the first place, two tautomeric forms of 
the molecule are possible by the addition of a 
proton to the xanthene nucleus, v7z. : 


HO O 0 oO O O 
~~ (NA : MA JN / 
= | H 

\4ZFH 0 YY 
c% | 7% 


| 
i ™ PS 
No 


oe) 


O— ZX 


a 


For each of these there can be written a corre- 
sponding structure, of higher energy: 


HOt O oo oe O OHt 
COO" CO 
\/ \ O 4\4\¥ O 
Ff | A: 
Cc Cc 
i ™ oN 
O O 
ai Bi 


The energy levels of the lower structures a and 8 
will not differ considerably. This will be true 
also with respect to the higher (dipole) structure 
a; and fi, so that we should expect an absorption 


TABLE I. Absorption bands of fluoresceins. 


























Am At (u.v.) 

Dye in mp - in mu R=\m/M 
Fluorescein 491 318 1.55 
Dichlorfluorescein 502 334 1.51 
Dibromfluorescein 504 328 1.54 
Diiodofluorescein 507 335 1.51 
Tetrabromfluorescein 518 
letraiodofluorescein 526 356 1.48 

TABLE II. Dissociation constants. 

Compound pK value 
Fluorescein pPKi=3-4.0 pK2=5.9 
Dichlorfluorescein (pK =4.1) 
letrachlortetraiodofluorescein (pK =4.3) 
Tetraiodofluorescein (pK =4.55) 





> A. Sklar, Rev. Mod. Phys. 14, 232 (1942). 
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Fic. 2. Spectral absorption of sodium erythrosinate. 
- = 1073 mole/I; 0.1 mm cell: X =10~4 mole/I; 1.05 mm cell. 
© =10-5 mole/1l; 10 mm cell. A = 10-6 mole/1; 100 mm cell. 
Only limiting curves showing dilution effect from hydrolysis. 


band from the transitions a—a,; and 6-8; fairly 
sharply defined, at least of the same order as 
that of the erythrosinate ion. However, in 
respect of the position of the band, i.e., the 
magnitude of hv for a—a:, B—6;, the energy 
levels of the a,8; configurations would appear at 
first sight to underwrite a band at higher 
frequency or shorter wave-length than that for 
the erythrosinate ion, where the two structures 
are of nearly equal energy. Somewhat surpris- 
ingly, the position of the band for the erythrosin 
acid (prepared by precipitation from aqueous 
sodium erythrosinate) appears to be in the same 
spectral region as that of the erythrosinate ion, 
at least in alcoholic solution, and indeed (in 
ethanol) at a longer wave-length, 545 my, than 
the ion (sodium erythrosinate) in water at 
527 mu. 

A somewhat remarkable result was obtained 
on comparing the absorption of sodium erythro- 
sinate in ethanol (nearly anhydrous) with 
erythrosin acid in ethanol saturated with NaOH. 
It was expected that the absorption curves would 
be identical, but as may be seen (Fig. 3) the 
sodium erythrosinate in ethanol gives a maxi- 
mum at longer wave-length than the erythrosin 
acid in ethanol saturated with NaOH, and 
practically identical with that of erythrosin acid 
in pure ethanol: 


Sodium erythrosinate in ethanol Am 544, 
Erythrosin acid in ethanol Am 545, 
Erythrosin acid in ethanol saturated 

with NaOH An 533. 
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Fic. 3. Spectral absorption of erythrosin. 1=Erythrosin 
acid in ethanol saturated with sodium hydroxide, 2 = Sodium 
erythrosinate in ethanol, 3=Erythrosin acid in ethanol. 


The behavior suggests that sodium erythrosinate 
in ethanol is converted almost entirely to 
erythrosin acid 


Na+ E=+2 EtOH—-EH2+2 EtONa. 


With increasing addition of water the absorption 
reverts to that given by the erythrosinate ion in 
water (see Table ITI). 

The longer wave-length of the erythrosin acid 
in ethanol (also in butanol) than of the erythro- 
sinate ion in water is further understandable by 
the following considerations. A greater influence 
of polarity of the solvent on the absorption 
would be expected for the erythrosin acid than 
for the erythrosinate ion—on the basis of the 
structures suggested—because a neutral struc- 
ture is supposed to be in resonance (to some 
extent) with a bipolar structure aa, and the 
latter should be favored by more polar solvents.* 
On the other hand, the structures of lower 
energy (a@ and 8) might tend to lose their hydro- 
gen-bonding stability in associating polar sol- 
vents, such as alcohols, thereby decreasing the 
energy difference between a and a, 8 and (Ay. 
The peculiar oscillation of the erythrosin acid 
absorption with solvent and pH may be caused 
by these factors, affecting both the tautomerism 
and the mesomerism. 

The low solubility of erythrosin acid in water 
is consistent with a high degree of intermolecular 


* The alcohols act with merocyanines as super-polar 
solvents. See reference 4. 
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hydrogen-bridging, and the extended absorption 
spectrum of the solid erythrosin acid (as sublimed 
film) seems in agreement with this (Cf. Fig. 4). 
At the same time, in the solid the maximum 
absorption appears to have shifted to the aux- 
iliary 6-band, a ‘‘dimerization”’ effect, which does 
not appear to occur to any great extent in the 
aqueous solutions of sodium erythrosinate.‘ Per- 
haps a closer packing side by side of erythrosin 
structures is allowed in the solid (acid) than is 
possible in the solution. 

According to the structures aa), B>6; sug- 
gested for the erythrosin acid, this could take on 
a third proton, giving the resonating structures 


HO* O O HO O OHt 
IVAN raf te \ (Np 
: — 
COD O COD O 
it it 
4 N\| NI 
O O 
(E~+3 H*) 
YI y2 

i and y2 shown above. Of these, y1 would tend 
to be somewhat more stable because of the 
possibility of double hydrogen-bonding as indi- 
cated. Thus, there would be de jure a crossing 
over of hydrogen-bonding stabilization from the 
6 to the a side, but not actually, because y; and 
Y2 are mesomers, not tautomers. The spectrum 
of erythrosin acid dissolved in sulfuric acid and 
alcohol, which is shown in Fig. 5, seems to agree 
with this. There is apparently a residue of the 
two-proton addition structure, indicated by the 
band at 540 mu, but in addition there appears a 
diffuse band having its principal maximum at 
487 my; the addition of a third proton has been 
markedly hypsochromic. 

There are two other points concerning the 
relations of the erythrosinate ion and erythrosin 
acid. The first is that the alkaline erythrosinates, 


TABLE III. Wave-length of maximum absorption of sodium 
erythrosinate in water-ethanol mixtures. 








Percent water Amax in mp 





0 540-544 
10 538 
50 533 
99 528 

100 527 
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e.g., the sodium salt, are somewhat hydrolyzed 
in water. This results in some deviation from 
Beer’s law on dilution, if not kept sufficiently 
alkaline, but it is not serious, and can be cor- 
rected. The other point is that in the exposition 
of structures for the acid form we have not taken 
account of deviations from coplanarity of the 
—COOH group with the xanthene residue. This 
is discussed fully in dealing with silver erythro- 
sinate, when it assumes a greater importance. 


2. THE ADSORPTION OF FLUORESCEIN DYES 
TO THE SILVER HALIDES 


Very extensive and solid foundations for our 
understanding of this field were laid by K. 
Fajans and his collaborators, in work extending 
over a period of years and from which they 
developed the important analytical technic of 
the “‘adsorption indicators.”’ In general, we shall 
not need to cite specific publications, since we 
can refer to two summary articles covering the 
field, by K. Fajans and T. Erdey-Griz up to 
1931° and by K. Fajans up to 1937.7 


The Initial Mechanism of Adsorption 


Various “‘adsorption rules’ are discussed by 
Fajans: There has been considerable disagree- 
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Fic. 4. Spectral absorption of sublimed erythrosin acid. 


°K, Fajans and T. Erdey-Griz, Zeits. f. physik. Chemie 
A158, 97 (1931). 
7K. Fajans, Die chemische Analyse 33, 161-222 (1935). 
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Fic. 5. Spectral absorption of erythrosin acid in ethanol 
acidified with 1 cc sulfuric acid per 100 cc ethanol. 


ment as to the relative preponderance in the 
initial adsorption act (a) of simple electrostatic 
attraction (of negative fluoresceinate ions to 
positive silver ions); (b) of exchange adsorption 
or displacement, e.g., of Br~ ions by fluores- 
ceinate; (c) of syncrystallization with the ad- 
sorbent ; (d) of the “‘solubility’’ of the adsorbate. 
As regards the adsorption of erythrosin to the 
silver halides, we may say that it is undoubtedly 
initiated by the electrostatic attraction of (posi- 
tive) Agt ions for the negative E= ions of 
(sodium) erythrosinate. However, the degree of 
adsorption (e.g., as measured for a_ specific 
surface by x/m—moles dye adsorbed per gram 
adsorbent) is not exclusively and definitively 
determined by the pAg (or pX) of the solution 
in the sense that the half-degree of adsorption 
coincides with the isoelectric point of the silver 
halide—as indicated in the theoretical curves of 
Fig. 6. These curves are constructed (a) for 
silver bromide, for which the equivalence point 
and the isoelectric points nearly coincide—as 
they do for silver chloride. The deviations of 
silver iodide will not be considered at present in 
this connection. (b) For silver erythrosinate, 
isoelectric point at pAg 5.12; (c) experimental 
results. 

Qualitatively, Fajans and his collaborators’ 
and Sheppard, Lambert, and Keenan*® have 
shown that the pAg determines the adsorption 
in the sense of the electrostatic theory. But it is 
believed that quantitatively there are deviations 
from the expected adsorption densities because of 
(a) the bulk factor of the dye ion. For example, 


8S. E. Sheppard, R. H. Lambert, and R. L. Keenan, 
J. Phys. Chem. 36, 174 (1932). 
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Fic. 6. Effect of pAg on adsorption of erythrosin. 
(a) Theoretical curve through J, AgBr, (b) Theoretical 
curve through 7, Ag2E, (c) Observed values for erythrosin 
adsorbed to AgBr. 


as discussed later, according to their orientation, 
a dye ion may be equivalent spatially to from 
three to ten or more of surface silver or halide 
ions. (8) Reversibility and solubility. Thus, the 
electrostatic hypothesis actually assumes a re- 
versible equilibrium in which the equilibrium 
constant depends upon the solubilities (or 
solubility products) of the adsorbent and ad- 
sorbate—and that these are independent of the 
fact of adsorption itself. Now the adsorption 
may not be reversible, because of orientation of 
the dye ion to present a non-lyophile aspect.‘ 
Again, the solubility of the adsorbed layer may 
be less than that of the substance in the mass, 
also probably the consequence of orientation. 
It is suspected that this is the case with silver 
erythrosinate. 

The existence of ‘‘exchange’”’ adsorption (8) is 
really tantamount to deviation from equilibrium 
determined by simple electrostatic attractions 
and the solubility product. From experiments 
on the adsorption of radioactive bromide ions, 
Fajans’ concluded that bromide ions are replaced 
by erythrosin, and are found in the solution, 
indicating that the solubility of ‘‘adsorbed”’ 
silver erythrosinate became equal or less than 
that of silver bromide. 

Something of the sort seems indicated by our 
experimental results in certain optical sensitizing 
experiments. Thus, a silver bromide emulsion at 
a pAg of 8.3 was found to show quite strong 
adsorption and optical sensitizing. The afore- 
mentioned equilibrium pAg measurements show 
that bromide ions must have been displaced. 
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It is true that the solubility product of silver 
erythrosinate seems to vary with the particle 
size (see later), but the lowest value, i.e., 
1X10-* mole/liter, is still much higher than that 
of silver bromide. At pAg~8, from the solubility 
products, silver erythrosinate should not be able 
to form on silver bromide. The explanation .of 
these anomalies may lie in the structure of solid 
silver erythrosinate, and in particular, of an 
adsorbed layer of this. 

Syncrystallization, in the sense of full iso- 
morphism, is excluded with dye ions like erythro- 
sin; from their magnitude and structure these 
ions are incapable of building in to the silver hal- 
ide lattice. 

The solubility factor has required a rather 
full study of the formation and properties of 
silver erythrosinate. According to a reversible 
electrostatic adsorption, at sufficiently low pAg 
values, the maximum adsorption densities of the 
fluoresceinates would be in the inverse order to 
the solubilities of their silver salts. Qualitatively, 
the data are in agreement with this. Thus, 
fluorescein and dichlorfluorescein do not show 
appreciable adsorption or sensitizing. However, 
in accordance with a reversible adsorption to a 
monomolecular layer, the magnitude of the 
adsorption approaches a limit with increasing 
concentration of erythrosin. 

The mathematical expression for this would 
be Langmuir’s equation for reversible formation 
of a monolayer, e.g., in simplified form: 


x/m=a=aBe,/(1+8c,), 


where x/m=a is the adsorption density, ¢, is 
the residual (equilibrium) concentration of the 
solution, and a@ and £6 are constants. Fajans’ 
states that the equation holds for erythrosin 
and eosin on silver iodide over the whole concen- 
tration range. Our results with erythrosin on 
silver bromide are not representable by this 
simple form of the Langmuir equation, but show 
approach to accord as high concentrations and 
adsorption densities are reached. 

The degree of fulfillment of the simple 
Langmuir equation is most readily exhibited by 
a plot of a against a/c,, which should be linear. 
In Fig. 7 may be seen the deviation of our results, 
and their approach to agreement at high ad- 
sorption densities. 
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An explanation for this may be found on 
consideration of the physical meaning of Lang- 
muir’s equation, as he originally developed it. 
Writing the equation in the form 


q=TH/(1+pa/Vi), 


where »=number of molecules (ions) adsorbed 
per unit area, r=average life of adsorbed ion 
(time of relaxation), .=g moles impinging per 
cm? per sec., a=fraction adhering per cm?, 
V,i=rate of evaporation or solution. It will be 
seen on comparison with the preceding empirical 
form that @ is proportional to 7, that c¢, is 
proportional to yw, and that simple relations 
obtain between the constants. It is now suggested 
that the discrepancy may be explained by sup- 
posing that 7, the average life of an adsorbed 
ion, is not constant and independent of the 
degree of adsorption in the case of erythrosin. 
We may suppose that 7 is affected by an ‘‘ac- 
commodation coefficient,’’ such that as the 
adsorption density increases, 7 decreases, ap- 
proaching a small constant value. This would 
be caused by the increasing difficulty of later 
adsorbed erythrosinate ions to accommodate 
themselves to the sparser sites, i.e., in attaching 
silver ions of the silver halide lattice to them- 
selves, and building in to the already formed 
monolayer of silver erythrosinate. This would 
give a curve deviating from the simple Langmuir 
function in the manner shown in Fig. 7. 

The determination of the limiting adsorption 
density (saturation) brings us to the effective 
area occupied by the dye ion and therewith to 
consideration of its orientation on adsorption. 
Because the solubility of silver erythrosinate is 
lowest for the dyes considered—and indeed (see 
later) appears to be effectively still lower for the 
adsorbed state—only this case need be con- 
sidered. The theoretical limit values for mono- 
layers, first calculated by M. L. Huggins® gave 
the following : 

Vertical 


Dye End-on Edge-on Flat Experimental 
Erythrosin 6 6 10 5.25 to 6 


The area covered by Evert in projection is about 


* Cf. S. E. Sheppard, R. H. Lambert, and R. D. Walker, 
J. Chem. Phys. 7, 265 (1939). 
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76A? and, as shown in Fig. 8, for Ena about 
144A”. It is evident from these that Ena: must 


cover >9 Agt ions, while Evert26 Ag* ions. 

The data of O. J. Walker!® gave E~23 Agt 
ions, which is definitely too high for any mono- 
layer. There seems to be no evidence for bi- 
molecular layers with erythrosin; it is probable 
that his determination of specific surface was 
considerably in error, since no size-frequency 
determinations were made of the microcrystals, 
but only random counts of linear measurements. 
Our own data for “‘saturation’’—which was not 
actually reached, because of independent forma- 
tion of silver erythrosinate—indicate a value of 
ca. 1 E= to 5~6 Ag*, using an area of 75A? for 
that occupied by erythrosin adsorbed edge-on 
or end-on. 

To secure more definitive ideas as to the 
nature of the orientation, it is necessary to 
consider the constitution of solid silver erythro- 
sinate, because the spectroscopic and sensitizing 
data show this to be fundamentally identical 
with erythrosin adsorbed to silver halide. 





afc 


Fic. 7. Adsorption of erythrosin to silver bromide. Curve 
A =Langmuir equation in linear form, Curve B=Com- 
posite curve of experimental data. 


10Q. J. Walker, Diss. Miinchen, 1926-27. 
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3. CONSTITUTION OF ADSORBED LAYER AND 
OF SILVER ERYTHROSINATE 


Upon adsorption of the fluorescein dyes to 
silver halides, there is a notable deepening of the 
color, compared with that of the ion in water. 
This deepening of the color is the basis of 
Fajan’s “adsorption indicator’? method and re- 
sults from a displacement of the absorption band 
to longer wave-lengths. The displacement from 
that of the (hydrated) erythrosinate ion is quite 
considerable, amounting to some 430A, or, as 
hv;—hv2, about=0.17 ev. Not only is that 
absorption band hy, responsible for the ‘‘ad- 
sorption indicator” function of eosin and erythro- 
sin, but it is also that effective in optical sensi- 
tizing."" As we have shown previously, it also 
corresponds to the absorption spectrum of silver 
erythrosinate. It has appeared desirable, there- 
fore, to investigate more fully the formation, 
properties, and structure of silver erythrosinate. 
This can be formed as a finely divided (micro- 
crystalline) precipitate of quite low solubility 
by mixing aqueous solutions of silver nitrate 
and sodium erythrosinate, followed by thorough 
washing and drying. O. J. Walker!® found that 
the wet precipitate tended to disperse colloidally, 
and only became microcrystalline after evapo- 
rating to dryness. He measured the solubility in 
water (probably near 20°C) by shaking for about 
24 hours with water, average value 1.1010-® 
mole/liter. Kieser! had previously obtained a 
value of 0.5X10-* and Pohl* 5.8X10-* mole/ 
liter. 

The formation and solubility product were 
investigated by electrometric titration of sodium 
erythrosinate with silver nitrate, accompanied 
by spectrophotometric measurements of the 
solutions and hydrosols. At the outset no special 
precautions were taken to exclude either light or 
carbon dioxide. Later, such precautions were 
observed, although any photolysis appears to be 
extremely slight or negligible. Since silver car- 
bonate has a solubility of about 3.7x10- 
mole/liter (25°C), it was possible that it might 
affect the titration, but later this appeared 
improbable. 


1S, E. Sheppard, R. H. Lambert, and R. D. Walker, 
J. Chem. Phys. 9, 96 (1941). 

12K. Kieser (O. J. Walker, Diss. Miinchen, 1926-27). 

* Cited by O. J. Walker. 
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Fic. 8. Erythrosinate ion superimposed flat on a 111 plane 
of a silver lattice. 


Forward and reverse titrations were made of 
one millimole of one component against number 
of milliliters of 0.02 N solution of the other. 
The neutrality point appeared for both curves at 
two atoms of silver per molecule of sodium 
erythrosinate, the two curves being nearly 
symmetrically distributed about this point. This 
corresponds to the composition AgE; it would 
follow that the dissociation constants cf the two 
Agt ions must be very nearly equal in magnitude, 
as was found for the pK values. This neutrality 
point may be regarded as indicating the silver 
ion concentration in equilibrium with solid silver 
erythrosinate, Ag2E. Assuming two silver ions 
(in solution), i.e., 

[AgeE J=2 Agt+E= 
(solid) 

the value pAg=4.30 should give twice the 
solubility, or S=2.5X10-5. However, the in- 
flection point was obtained from the differential 
curves at pAg=5.04, which, on the above 
assumptions, gives S=4.55X10-* mole/liter at 
25°C. 


Measurement of Solubility of Silver Erythro- 
sinate by Dilution of Suspensions 


Suspensions of silver erythrosinate were made, 
which differed in respect of excess silver or 
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erythrosinate ion, and the pAg was measured. 
These suspensions were then successively diluted 
with water, and pAg was measured at each stage. 
Curves were plotted for log (Ag2E) against pAg 
observed. 

Solutions of pAg>5.15 decreased in pAg 
through a minimum then increased with further 
dilution. Those of pAg<5.12 increased with 
similar dilutions to pass through an inflection 
point at pAg 5.12. This value would imply a 
solubility of 3.8X10-® mole/liter. Ideally, the 
point at which all curves intersect should give a 
value on the dilution axis corresponding, as log S, 
to the solubility. This value appears to be at 
4.6-4.8, but is not sufficiently well defined to be 
considered with confidence. 


Spectrophotometry of Suspensions and Solutions 


Successive dilutions were made of a suspension 
3.24X10-* molar silver ‘erythrosinate (AgeE). 
Absorption curves were measured, and also pAg 
and pH measurements. The data are shown 
graphically in Fig. 9. During the time of exami- 
nation, the ‘‘solutions’’ all appeared clear, so 
that fairly accurate absorption coefficients could 
be measured. However, on standing 5 to 6 days, 
all specimens, up to and including that of 
4X10-* mole/liter of silver erythrosinate, had 
precipitated, so that evidently they had not 
reached equilibrium at the time of measurement. 

They show, however, the band at ca. 568 mu 
of a hydrosol of silver erythrosinate giving place 
on dilution to the band at 528 mu of sodium 
erythrosinate in solution. The data are also given 


toxi0°5 






09 
os 
o7 
06 
os 
04 
o3 


EXTINCTION COEFFICIENT 


a2 


0. 


° 
400 80 90 500 10 20 30 40 50 6O 70 80 90 G00 10 20 
WAVE-LENGTH IN MILLIMICRONS 


Fic. 9. Effect of pAg on spectral absorption of erythrosin 
X =3.24 10-4 N AgNOs, © =3.24X10—, O =1.62 10-5. 
O=8.1X10-§, A=3.24X10-§, -=1.62K10-7, ~=6.5 
x 1077. 
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in Table IV to bring out the behavior in respect 
of pAg. Other similar sets show that the pH 
values are fairly constant. 

Another similar experiment was made with a 
suspension or hydrosol initially 1.610-5 mole/ 
liter of silver erythrosinate; the band at 568 mu 
became indiscernible between pAg 5.33 and 
pAg 5.55. It is evident that silver erythrosinate 
has a strong tendency to occur and persist in 
colloidal solution. It is probable that this affects 
the solubility, making it greater for small parti- 
cles. This was corroborated by observations on 
sols freshly prepared, with varying ratio of silver 
ion and erythrosin, compared with identical 
systems kept agitated for two weeks by bubbling 
nitrogen through them, when spectrophotometric 
measurements were made on the filtrates (Table 
V). Values obtained directly from the inflection 
point on titration, viz., 4.5; 10-*® mole/liter, and 
by dilution of less ‘‘aged”’ sols or precipitates (e.g., 
4.5,X10-*, 3.8,10-*) fall intermediately be- 
tween these extremes. It appears definite that on 
“aging” a precipitate with agitation in water, 
the solubility falls, and may well approach that 
obtained by O. J. Walker for a well-washed 
precipitate evaporated to dryness, viz., 1.0 10~® 
mole/liter. It is possible that the solubility of 
the adsorbed layer is even lower, at any rate 
the behavior of erythrosin adsorbed to silver 
bromide indicates this. Taking the solubility of 
silver bromide as 0.85X10-* mole/liter, this 
would not differ greatly from that of O. J. 
Walker’s value for silver erythrosinate. Fajans’ 
observed that radioactive bromide ions of— or on 
—silver bromide could be displaced by erythrosin 
into the solution, which would point to the 
solubility of adsorbed silver erythrosinate being 
=Sa,pr; and as noted we have observed strong 
adsorption of erythrosin at pAg 8.3. It is possible 
that we encounter here an effect of orientation of 
the adsorbed erythrosin. As noted in Section 2, 
the orientation appears to be either edge-on or 
end-on, and some consequences of this will be 
discussed now in connection with the constitution 
of solid and of adsorbed silver erythrosinate. 

We have seen that no separable and soluble 
ion of the composition AgE~ exists. The band 
at 567-570 my is only produced by [AgeE ] 
adsorbed or solid and is an example of what 
have been termed “‘cooperative’”’ spectra. The 
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large displacement from that of sodium erythro- 
sinate was attributed by Fajans to ‘‘deforma- 
tion” of the erythrosinate anion by the large 
silver cation. If the normal structure of the 
erythrosin anion is considered (cf. Section 1), it 
is reasonable to suppose that cations such as 
Nat would tend to place themselves symmetri- 
cally with respect to the resonating oxygen atoms: 


hy P 
7 Na? (1) 
*\ 
VKH, O 
O- O 
" 
Oo 
N\ 
6A Nat (2) 


There does not appear to be any reason why 
substitution of an Agt ion for Nat (2) [still less 
for Nat (1) ] should change the resonance system 
so as to displace the absorption by ca. 0.2 ev. 
The Agt (2) ion will be at too considerable 
distance—of the order of 5 to 6A—to have any 
considerable polarizing effects on either O-atom. 
And, otherwise, it appears that the resonance 
system of the molecule-ion is shielded from any 
influence of a (coplanar) Ag* ion. 

Such an arrangement would presume an edge- 
on orientation of the erythrosin ion relative to 
an octahedral plane of the silver halide. We must 
examine the possibilities of an end-on orientation, 
since the adsorption-density data do not allow 
any preference. Another ‘“‘coplanar” arrangement 
for a silver erythrosinate molecule now suggests 
itself based on a possible assimilation to erythro- 
sin acid (cf. Section 1). Consider the following 


TABLE IV. Displacement of absorption. 











Concen- 

tration Log 
moles/liter concentration 

of Ag* calculated pAg ©5097 €567 
1.59 10-5 —4.80 4.80 0.264105 0.563105 
1.181075 —4,.927 5.04 0.268 0.543 
8.0 «10-6 —5.10 5.09 0.323 0.438 
6.3 1076 —5.20 5.16 0.337 0.337 
3.16 X 1076 —5.50 5.33 0.463 0.115 
1.59 10-6 — 5.80 5.55 0.474 0.060 
8.0 x10-7 — 6.10 5.78 0.485 0.050 
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graphic formula: 
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If we conceive that the ‘‘phenolic’’ Agt ion 
might—as with the H-band of erythrosin acid— 
extend the resonance system through the “car- 
boxyl” Agt ion—both silver ions being now 
“chelated’”’ between oxygen atoms in resonance 
circuits—we have a feasible explanation of the 
displacement to longer waves of the absorption 
—by extension of the resonance system as in the 
case of erythrosin acid. Such a constitution 
agrees also with the factual equality of the two 
Ag* ions in respect of dissociation, as shown by 
the titration curves. It must be admitted, 
however, that the completely coplanar configu- 
ration as shown above is very unlikely, if not 
impossible. The carbon atom of the carboxyl 
group could hardiy come as near to the 4 —CH, 
while with normal 125-degree bond angles the 
upper O-atom of the carboxyl could not be in 
the position assigned, but must, through rotation 
of the dependent phenyl, be in a plane consider- 
ably above or below that of the xanthene 
residue. This would be enforced the more because 
of the bulk of the iodine atom (—C—I=2.1A) 


TABLE V. Solubility. 











Age Solubility product Solubility 
Immediate 7.6X 10716 5.74X 10-6 mole/liter 
Two weeks 

with agitation 0.7 X 10716 2.5 X10-§ mole/liter 








attached to the —C atom adjacent to the 
“phenolic” group. These conditions obtain, how- 
ever, in much the same degree for erythrosin 
acid. Rotation of the secondary or accessory 
resonance system—involving the carboxyl group 
and the two Agt ions—to a plane making an 
angle of ca. 50 degrees with the xanthene residue 
might allow the distance between the two Agt 
ions to be increased to a reasonable extent. In 
a coplanar configuration they would be separated 
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only—between centers—by some 2.2A, which is 
too close. And this also disagrees with their 
probable interval as dictated by the silver halide 
lattice for adsorbed erythrosin. 

Admitting this non-coplanarity of the aux- 
iliary resonance system (involving the two Agt 
ions) there still remains another difficulty. With 
this end-on orientation of the silver erythrosin 
molecule, the principal transition moment must 
still be in the xanthene residue—as between the 
two quasi-phenolate oxygen atoms. But this 
means—for erythrosin adsorbed to silver halide 
crystals—that the vibrations absorbed must be 
polarized in directions more or less perpendicular 
to the plane of the crystal, that is, parallel to the 
axis of propagation of the light. This would imply 
longitudinal vibrations in the light ray, which 
are not admitted as possible. 

Whether the part of the suggested resonance 
system which would be in the plane of the crystal 
—and, therefore, transverse to the Poynting 
vector of the incident light ray—could be re- 
garded as taking care of this we cannot say. 
The configuration offers one possibility, we will 
consider another. It has been noted that the 
absorption band of silver erythrosinate is that 
of the solid phase. This suggests the possibility 
that in the solid, where a new order of aggrega- 
tion occurs, the absorption involves this order, 
and, therefore, more than one molecule. A 
possible structure of aggregation is shown in 
Fig. 10. In this arrangement, two erythrosinate 
ions are supposed in parallel, and separated some 
5-6A, both because of the non-coplanarity of the 
dependent phenyl with its carboxylate group, 
and also because of the diameters of the —CI 
groups. The new unit, composed then of two 
erythrosinate and two silver ions* is capable of 
resonance between two configurations of identi- 
cal energy, in which the interstitial silver ions 
act as alternate acceptors and donors of electrons. 
The absorption act would involve no permanent 
reduction of silver ion to metallic silver, which 
corresponds with observations on the optical 
stability of silver erythrosinate and with the 
observations on the photolysis of silver bromide 
sensitized with erythrosin. In the presence of 
halogen acceptor, the dye is not destroyed by 


* Neglecting the silver ions attached to carboxyls. 





the absorption of light which produces metallic 
silver from the silver halide." 

As a linear mesophase, this structure could be 
continued indefinitely,‘ and such a structure may 
occur in the monolayers of adsorbed erythrosin. 
But it could hardly be present in solid crystals of 
silver erythrosinate, because it involves a ratio 
of ‘‘phenolic’”’ silver ions to erythrosinate ions 
varying continuously with the length of the 
chain, and this is not compatible with the 
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Fic. 10. Suggested structure of adsorbed silver 
erythrosinate. 


stoichiometric composition of the crystal. How- 
ever, this does not affect the identity of the 
absorptions of the adsorbed layer and of the 
solid crystal, because the new absorption band 
is attributed solely, as regards wave-length, to 
the structure of the optical doublet or dimer. 
It is evident that this conception of the structure 
of the absorber for the new band of silver 
erythrosinate, and of erythrosin adsorbed to 
silver halides, must remain in the realm of 
hypothesis until accurate x-ray diffraction data 
for the crystal, and electron diffraction data for 
the adsorbed layer, become available. 

Another possibility for a participation of silver 
atom ions in a new resonance system was 
suggested by M. L. Huggins, in view of the 
arrangement shown in Fig. 11. This also involves 
charge transfer through silver atom ions located, 
however, between O--atoms of two successively 
aligned erythrosinate ions. A possible objection 
to this is that it might involve dependence of 
the absorption maximum upon the number of 
molecules aligned, i.e., on the length of the 
chain, and there is no evidence of this. 

The arrangement previously discussed has 
evident affinities with one suggested for the 
mesophase of cyanine dyes‘ and would refer light 
absorption and optical sensitizing by eosin and 
erythrosin and the like to similar polymolecular 
structures. It is possible that the photochemical 
(quantum) yields in the photolysis of sensitized 
silver halides may give another datum deciding 
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this point. Should similarity of sensitizing by the 
fluorescein type of dye with that by the meso- 
phases of cyanine dyes be established, we could 
discriminate three main classes of optical sensi- 
tizing by dyes: 


(1) By neutral (non-ionized) dye molecules: 
e.g., by merocyanines, by 1,4-diamino- 
anthraquinone (and other anthraquinone 
derivatives) by porphyrins, as chlorophyll, 
phthalocyanine, haemin. 

(2) By (basic) cyanine dye molecules, and by 
thiazines. 

(3) By the mesophases of cyanine dyes, and 
by the fluorescein dyes. 


Such a reduction to specific types of absorbing 
units could give a greater precision to our notions 
of the mechanism of optical sensitizing of the 
silver halides. 


4. OPTICAL STABILITY AND SENSITIZING 


Previous experiments were made with dry 
silver erythrosinate powder, and with silver 
erythrosinate suspended in gelatin (de-ashed and 
thoroughly freed from halide ions), washed with 
distilled water, and dried. Only with very intense 
illumination (restricted as well as possible with 
filters to its own absorption band) could any 
reduction to silver be observed, and this was 
attributed to the heating effect. Normal up to 
large exposures, with the gelatin-emulsified silver 
erythrosinate, failed to give any latent image 
developable with regular chemical developers, 
although prefixation physical development was 
possible. 

In connection with the experiments on the 
constitution of silver erythrosinate (Section 3), 
some decomposition of silver erythrosinate sus- 
pension was observed on long exposure to day- 
light, darkening and separation of silver occur- 
ring. These were only over periods of time— 
days—entirely incommensurate with normal 
photochemical sensitivities, and occurred only 
with excess of silver ion. It would be attributed 
to optical sensitizing of adsorbed silver salts. 

Silver erythrosinate must be regarded as 
properly an optically stable absorber in the sense 
that it is per se able to dissipate the energy taken 
up by absorption of light as heat or fluorescence. 
This is in good accord with its behavior as an 
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Ag*ion plane of AgBr crystals 
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optical sensitizer for the silver halides. It was 
with erythrosin-sensitized gelatino-silver bromide 
emulsion that Lescynski' first made the basically 
important observation that one sensitizing dye 
molecule could produce in photolysis a large 
number of silver atoms." 

The persistence of the sensitizing molecule 
(for silver halide sensitized with silver erythro- 
sinate) only occurs in the presence of a halogen 
acceptor. For these reasons we consider that the 
mechanism of optical sensitizing and regenera- 
tion of (silver) erythrosinate proposed by S. J. 
and G. L. Natanson™ is inacceptable. They 
suggest the obviously tempting view that a 
negative charge (electron) is irreversibly trans- 
ferred from the erythrosin part of the complex 
to the silver ion part. The chemically changed 
erythrosin* is supposed to recover an electron, 

18 W. Lescynski, Zeits. f. wiss. Phot. 24, 275 (1926). 

14 Cf. S. E. Sheppard, R. H. Lambert, and R. D. Walker, 
J. Chem. Phys. 7, 426 (1939). 

15S. J. and G. L. Natanson, J. Phys. Chem. U.S.S.R. 14, 
278, No. 2 (1940). 


* The primary products of such electron loss would be 
free radicals. 
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either from bromide ion, or from “reducers” 
such as sodium nitrite—or from gelatin, which 
is hardly a reducer (of silver ions)—but is a 
bromine acceptor. 

Bromine, as gas, is rather slowly taken up by 
erythrosin, with some change of color. There 
appears to be a twofold reaction, in that an 
irreversible attack on the erythrosin occurs, but 
also some substitution of bromine for iodine in 
the molecule. 

In the presence of water, bromine forms 
hypobromous acid, according to the reaction 
Br2+H,O=HOBr+HBr (which is shifted to 
the right by hydroxy] ions, the reaction analyzing 
to Br2+2 OH-=OBr-+Br-+H,0). A protec- 
tive acceptor, in optical sensitizing, should, 
therefore, not only be a halogen acceptor, but a 
reducer reacting with the hypobromite, which 
can destroy the dye. There is evidence!® that 
gelatin, in the ordinary photolysis (not optically 
sensitized) of gelatino-silver bromide, reacts 
effectively with hypobromite to protect the 
latent image silver. 

It is probable that the erythrosin molecule 
can act as a temporary or intermediary bromine 
acceptor by bromine adsorption to the iodine 
atoms of the molecule. Given sufficient time, and 
nothing to react with the bromine, there is 
displacement of iodine by bromine—the cuckoo 
in the sparrow’s nest !—but this probably plays 
no part in optical sensitizing of silver halides. 
In this case the bromine may be temporarily 
held, only to react later with water and with 
halogen acceptors proper. 

The balance of the evidence is that, in optical 
sensitizing, the silver erythrosin complex on 
silver halide is not chemically affected by the 
light it absorbs, but transfers energy as an 
exciton to the silver halide, with which it is 
combined by “adsorption.” In this primary 
phase of sensitizing, the mechanism does not 
differ from that with basic cyanine dyes. These 
are positive ions, and the emission of an electron, 
as suggested by various authors!’ is much more 
difficult to conceive. 

16 Cf. S. E. Sheppard and C. L. Graham, J. Frank. Inst. 
230, 619 (1940). 

17 Cf. G. Kornfeld, J. Phys. Chem. 42, 795 (1938); also 


S. E. Sheppard, Atti. X. Congr. Internaz. Chim. Rome 1, 
234 (1938). 
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Discussion of Energetics 


The adsorption of erythrosin may be conceived 
as implying the participation of silver ions of the 
silver halide lattice in the resonance system of 
the erythrosinate. In this sense we may regard 
the resonance of the new system as equivalent 
to AgEAg* E- (disregarding for the moment 
the second electrovalence). We do not know 
the energies of the two states—of the homopolar 
and dipolar configurations, respectively, but it 
is legitimate to suppose that of the dipolar form 
is higher. Absorption of a quantum of light 
results in an excited state of AgtE~-. In the 
solid unadsorbed silver erythrosinate, the exciton 
is dissipated!* by degeneration to thermoelastic 
vibrations with relapse of the system into the 
AgE condition. When adsorbed to a silver halide, 
the exciton—polarized in the plane of the 
excited absorber—is transferred to the silver 
halide lattice and results in the neutralization 
of a particular pair of Ag+ X~- ions—probably in 
an ‘‘altered’’ region of the lattice. It is reasonable 
to assume that there is such a region immediately 
contiguous to the adsorbed dye. The dye has 
been stabilized by adsorption (as shown by its 
lower solubility than that of the ‘‘pure’’ precipi- 
tated material, vide infra) at the expense of the 
stability of the silver halide lattice. The syn- 
polarization of an affected Ag+ Br~ pair has been 
considered previously." Where a ‘‘sensitizing 
nucleus’”’ exists, the exciton, handed on by the 
resonance Agt X-AgX, finally provokes an 
irreversible separation to Br+(Agt+ @)—giving 
a quasi-free, or photoelectron. There appears 
nothing in this inconsistent with photoconduct- 
ance, nor with the Gurney-Mott theory of 
photolysis. Such a scheme might assume equiva- 
lent migrations of electrons and positrons (nega- 
tive potential holes), but effectively it is con- 
venient to assume that positive charges remain 
attached and move with silver atoms as silver 
ions—while the electrons move relatively freely. 

The question as to whether adsorbed dye and 
silver halide lattice have preponderantly elec- 
trons in common—as part of a new complex— 
or, have excitons in common, will be discussed 
later. 


18 J. Franck and E. Teller, J. Chem. Phys. 6, 861 (1938). 
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The Transmission Coefficient in the Theory of Absolute Reaction Ratest 
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The activated complex theory of reaction rates is discussed critically. The transmission 
coefficients for a number of idealized energy surfaces have been computed quantum mechani- 
cally. It appears that curvature of the reaction path may introduce activation energy in 
addition to that caused by constriction and elevation of the energy valley. Except at low 
temperatures, this additional activation energy is almost negligible. However, it may contribute 
to the separation of isotopes, where small differences in activation energy are important. 
For systems in thermal equilibrium at room temperatures and above, the average quantum 
mechanical transmission coefficient is found to differ insignificantly from that calculated by 


classical mechanics. 





N the theory of rate processes proposed by 
Eyring,! the process is represented in configu- 
ration space by the classical motion of a particle 
on a surface, the height of which is proportional 
at each point of the surface to the potential 
energy of the reacting system whose configura- 
tion is represented by that point. Such surfaces 
have the general shape of long, narrow valleys 
joined by ‘‘passes” through the regions of high 
energy. The strip of configuration space lying 
across the deepest of such saddles connecting two 
regions which correspond to the initial and final 
states of the system is called the activated state. 
Whenever the motion of the representative par- 
ticle on this surface can be treated classically, the 
rate at which systems move from initial to final 
configurations is given by 
P t d 
k=—-k (1) 
P;6 

and is a direct consequence of statistical me- 
chanics. In this equation, P; and P; are the prob- 
abilities of the initial and of the activated states, 
respectively ; 6 is the width of the transition state; 
v, the average velocity of the particles moving 
across 6; and x is the transmission coefficient, the 
fraction of particles successful in completing the 

journey across the activated state. 

{+ Taken in part from a thesis submitted by H. M. 
Hulburt to the Graduate School of the University of 
Wisconsin in partial fulfillment of the requirements for the 
degree of Doctor of Philosophy. — : 

Present address: Frick Chemical Laboratory, Princeton 
University, Princeton, New Jersey. 
** Present address: Geophysical Laboratory, Washing- 


ton, D. C. 
1H. Eyring, J. Chem. Phys. 3, 107 (1935). 
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Two basic assumptions are involved in this 
theory: 

1. Equilibrium is maintained between the ini- 
tial and activated states.throughout the course 
of the reaction. 

2. Classical mechanics is adequate to describe 
the motion of the reacting systems over the entire 
course of the reaction. Both of these assumptions 
have been the subject of much discussion.? The 
first assumption is justified a posteriori by the 
success of the theory, but it has been difficult to 
see in advance the conditions under which it 
holds. As to the second assumption, it is evident 
that classical mechanics is not applicable to all 
the motions of the system as it passes through 
the activated state. A strictly rigorous treatment 
involves replacing the classical particle by a 
quantum mechanical wave-packet subject to the 
potential field represented by the energy surface. 
It seems probable at first that there should be 
significant deviations from the behavior pre- 
dicted classically because of reflection, diffrac- 
tion, and tunnelling of the quantum mechanical 
wave packet. Such effects might be expected at 
barriers, constrictions, and bends in the valley in 
the energy surface in which the reaction path 
runs. However, it has previously been shown? 
that for most chemical reactions, the reflection at 
barriers and constrictions in a straight path is not 
great enough to modify the classical predictions 
significantly, except possibly for reactions at very 


low temperatures, or for reactions involving iso- 


2 (a) Summarized in Fowler and Guggenheim, Statistical 
Thermodynamics (Cambridge, 1939) p. 517; (b) J. O. 
Hirschfelder and E. Wigner, J. Chem. Phys. 7, 616 (1939). 
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tope separations. The effect of curvature in the 
reaction path has not been considered previously 
and will be discussed here. 


QUANTUM EFFECTS 


The magnitude of the quantum effects in the 
motion of the system in curved reaction paths 
can be estimated by calculating the transmission 
coefficient for a wave packet traveling in the 
given potential field. The mathematical difficulty 
of treating the problem restricts our attention 
to a few highly idealized potential surfaces, from 
which, however, we can discern the general fea- 
tures of the more complicated cases. So as to 
present the results as vividly as possible, most of 
the detailed mathematical development will be 
given in the appendix, to which reference may be 
had for the justification of the qualitative argu- 
ments presented here. 

It is important to distinguish between two 
types of reaction at a barrier or bend in the 
energy surface. If the electronic quantum num- 
bers of the reacting complex do not change 
throughout the course of the reaction, the reac- 
tion is called quasi-static. The activated state 
method is restricted to such reactions. There are 
then two possible types of reaction: 


1. that in which the vibrational quantum 
numbers of the reacting complex remain un- 
changed during the course of the reaction; 

2. that in which vibrational transitions take 
place. 


The conditions for good transmission are very 
different in the two cases. 


A. No Vibrational Transitions 


We shall consider first the case in which vibra- 
tional quantum numbers are unchanged by re- 
flection or transmission at an irregularity in the 
energy surface. When this is the case, the Eyring 
formulation of the theory of absolute reaction 
rates is easily justified for incident systems with 
sufficiently great translational energy. Mathe- 
matically, the condition for no change in quan- 
tum number is realized when the potential energy 
is such that the Schrédinger equation is separable 
in coordinates parallel and perpendicular to the 
reaction path. The wave functions for vibrational 


and translational motions aré then completely 
independent. Such separability is never complete 
in any actual case, but it will be approximately 
attainable if the vibrations perpendicular to the 
reaction path are rapid compared with the mo- 
tion parallel to the reaction path. The system 
then executes many vibrations before moving 
into a region where the vibration frequencies are 
much different. Hirschfelder and Wigner? have 
shown that, for a straight reaction path, if the 
change of vibration frequency, v, along the reac- 
tion path, s, satisfies the condition 
dv 


v—<y" 
ds 


where v is the translational velocity of the system, 
there will be no change in the vibrational quan- 
tum number. 

The transmission coefficient for the case of no 
change in vibrational quantum number where the 
reaction path is curved may be studied by con- 
sidering the potential energy surface to be the 
channel between two confocal hyperbolae. Such 
a channel becomes narrower as the system ap- 
proaches the greatest curvature and then widens 
out. The potential energy is zero inside the chan- 
nel and infinite outside. For this case the Schréd- 
inger equation is separable and can be solved ex- 
plicitly for the whole region. Hence, there is no 
possibility for a change in vibrational quantum 
number in rounding the bend. By introducing the 
elliptic coordinates, u and ¢ such that 


x=R cosh yu sin ¢, (2) 
y=R sinh p cos o 


we can transform the Schrédinger equation with- 
out potential energy into the form: 


OY Fy 82°m RE 
—] (cosh 2u+cos 2¢)y=0. (3) 








0 ue 0¢? h? 





The lines ¢=¢o are the hyperbolae (x?/sin? ¢o) 
— (y?/cos? ¢9) = R? whose asymptotes have the 
slope dy/dx=cot ¢. Hence, ¢o is the angle be- 
tween the asymptote and the y axis. Let the 
angle between the asymptotes of the two bound- 
ary hyperbolae be A. Then the boundary condi- 
tion is that y=0 on the two hyperbolae whose 
asymptotes are defined by @=¢9+A/2 and 
¢=¢)—A/2. Equation (3) is separable into the 
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two one-dimensional equations, 
@U 8nx°m 

ae he 
ab 81m 

dg ht 


where ¥(u, ¢) = U(u) -®(¢). Here \ is the separa- 
tion constant whose value is to be determined by 
the boundary conditions. 

Since the lines u= yo are ellipses, (x?/cosh? yo) 
+(y?/sinh? wo) = R? orthogonal and confocal to 
the hyperbolae, ¢=¢o, the factor ® gives the 
variation of the wave function y as the system 
moves along any ellipse, i.e., across the channel, 
while the factor U gives the variation as the sys- 
tem move along the channel, parallel to @o. It is 
evident that for a given value of » and ¢ there 
are two values of x and y, viz., (x, +y) and 
(x, —y). Hence, the wave function is symmetrical 
about the x axis. Its derivatives, however, are 
different for y positive and for y negative. In 
order to make them equal at y=0, it is necessary 
to add a component representing a reflected wave 
to that for the incident wave in the region y>0. 
Thus it is necessary to introduce a reflected wave 
to satisfy the continuity condition on the wave- 
function. Physically, however, this reflection is 
most easily interpreted in terms of a virtual 
barrier, which can be demonstrated without 
actually solving the equations (4) and (5) 
completely. 

We solve Eq. (5) for , determining the allowed 
values of \. When these values of \ are substi- 
tuted into Eq. (4) it becomes a one-dimensional 
equation for the translational motion of the sys- 
tem through the channel. For the lower energies 
and for values of A less than 30°, the equation (5) 
may be solved using an asymptotic expansion for 
® and the allowed values of \ can be developed 
in a power series in A. (See Appendix). We then 
find for \, up to terms in A?, 





E 
R(5 cosh 2u-n) U=0, (4) 





E 
R(> cos 26+n)o=0, (5) 





a a 
= —— cos 2¢0— 
8mR2A? 2 : 


where b(£) is a function of the energy given in 
the Appendix. 

Putting Eq. (6) into Eq. (4) gives a one- 
dimensional equation for the motion of the par- 
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ticle parallel to any hyperbola. This equation 
is not in the usual Schrédinger form, how- 
ever. It can be put into the usual form by the 
substitutions, 








" 
s=R{ (sinh? u+cos? ¢o) ‘du, (7) 
0 
x = U(sinh? u+cos? ¢o)?. (8) 
Equation (4) then becomes 
Seti Niews 
ds? ft? " — 
where 
h? 
Vs= 
822m R? 





[= 3 cosh? 2u+2 cosh 2u cos 2¢ 


10 
(cosh 2u-++cos 2¢0)* (19) 


and 


V h2n2 soce)|/ 
* L8mR2A? 


3(cosh 2u+cos 2¢9). (11) 





The coordinate s, defined by Eq. (7), is the dis- 
tance from the narrowest point of the channel 
measured along the hyperbola asymptotic to 
x=y tan ¢o. Thus, s is the coordinate of the reac- 
tion path lying midway between the walls of the 
channel. The substitution of Eq. (8) is made so as 
to bring the equation into the form (9), where x 
now describes the translational motion of a par- 
ticle in one dimension, moving in a virtual poten- 
tial field which is the sum of V, and V4. 

The separation of the virtual potential into 
two terms is justified by comparing V, with the 
energy, Evi», of the mth vibrational state for a 
particle confined to a box of width d given by 


A 
2=4R?(sinh? u+cos? ¢o) sin? > (12) 


We see then that V, is just equal to Ey for 
values of A small enough that terms in A? and 
higher powers can be neglected. Thus we find 
that though the height of the energy surface is 
perfectly constant, a system moving in a curved 
path encounters an effective potential energy 
barrier, which is partially due to the change in 
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Fic. 1. Elliptic channel. Classically, all particles whose 
trajectories lie in the angle a are transmitted. 


zero-point energy of vibration as the channel 
narrows (V,) and partially due to the curvature 
of the path (Vz). 

If the channel is elliptic rather than hyperbolic, 
the zero-point energy due to the constriction in 
the channel has a close connection with the class- 
ical conditions for transmission. In the classical 
case, the particle is reflected in straight lines from 
the walls of the channel. At any point, its tra- 
jectory is specified by the ratio of vibrational to 
translational energy, for we have (Eyi»/Etrans)* 
=v, /v,,;=tan 6 where v, is the component of ve- 
locity perpendicular to the reaction path and »,, 
is the component parallel to the reaction path. 
Inspection of Fig. 1 shows that all those particles 
whose trajectories are such that @<a are trans- 
mitted, whereas those whose trajectories have 
6>aare reflected. Here there is, even in the class- 
ical case, a minimum translational energy for 
each vibrational level which must be exceeded 
before transmission occurs. An exact calculation 
shows that precisely the same critical value of 
(Eviv/Etrans)! is predicted in the quantum me- 
chanical case. Thus, for an elliptic path, there 
is a classical analogy for the virtual zero-point 
energy due to a constriction in the channel. 

The term V4 merits closer consideration. It is 
to be noticed that at large distances, V4 corre- 
sponds to an attractive potential, but nearer the 
narrow part it becomes a repulsion. This peculiar 
behavior may be interpreted with the help of 
three analogies from other branches of physics. 
The presence of Vz arises mathematically from 
the ‘straightening out’’ of the hyperbolic reac- 
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tion path when the transformations (7) and (8) 
were made. This result is analogous to the rela- 
tivistic description of gravitational forces in 
terms of geodesics in curved space. Thus a par- 
ticle experiencing gravitational force when its 
motion is described in euclidean space appears to 
move on a geodesic in no field in the curved rela- 
tivistic space. 

A second analogy can be made with the clas- 
sical motion of a particle constrained to move ina 
prescribed region. Such a particle experiences an 
apparent ‘‘force’’ due to its constrained motion. 
For example, the motion of a particle moving in 
an elliptical orbit may be described by supposing 
an effective central force, directed toward one 
focus of the ellipse. Similarly, here, the walls of 
the channel restrict the particle to move in a 
definite region of space. This gives rise to a 
virtual field of force acting on the particle. 

A third analogy makes use of the wave proper- 
ties of matter. A parallel beam of light is made to 
converge to a point by reflecting it in a parabolic 
mirror. Similarly, any curved mirror tends to 
focus waves, concentrating intensity at some 
point, and diminishing it at others. This focusing 
effect bears a strong resemblance to the quantum 
effect observed here, for in regions of low virtual 
potential energy, one should expect to encounter 
a concentration of systems, due to the Maxwell- 
Boltzmann distribution. The presence of curved 
reflecting walls appears to focus the systems to 
some extent, concentrating them in the regions 
where V, is low and diminishing their numbers 
where V, is high. 

The magnitude of V, is indicated by the size 
of the factor h?/8x’mR?, where R is the focal 
length of the hyperbola ¢=¢o. For the reaction 
H+H2, R~1A and M=? in atomic weight units. 
For D+D2, R is unchanged, but M=4#. Hence 
Vz will be only one-half as great. The maximum 
value of V5 is 0.288 kcal./mole for the H+A, 
reaction where ¢9= 60°. This difference is of the 
right order of magnitude to account for the dis- 
crepancy between the observed and calculated 
rates.® 

We may conclude that when there is no change 
in vibrational quantum number, we should ex- 
pect fast particles to be transmitted more rapidly 


3L, Farkas and E. Wigner, Trans. Faraday Soc. 32, 1 
(1936). 
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a-7 Fic. 2. Classical transmission 
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than slow ones. This would appear experimen- 
tally as an additional activation energy beyond 
that necessary to surmount the energy barrier 
were the reaction path straight. 


B. Reflection with Change in Vibrational 
Quantum Number 


Thus far we have investigated only cases where 
the system is transmitted or reflected without 
change in vibrational quantum number. For a 
general potential surface this will be the case, as 
we saw, if the zero-point energy is a sufficiently 
slowly varying function of the reaction coordi- 
nate. This means that the actual slope of the 
reaction path must not be too great at any point, 
and that its curvature must be small. Even then, 
the curvature of the path will introduce an addi- 
tional activation energy for the reaction. In gen- 
eral, when these conditions are not satisfied, there 
will be reflection of the incident system and simul- 
taneous change of its vibrational quantum num- 
ber. Exceptional conditions prevail when the 
potential energy surface is such as to render 
the Schrédinger equation separable everywhere. 
Then there can be no change in quantum number, 
but reflection will still occur, due to the increase 
of zero-point energy as the channel is constricted 
and to the curvature of the reaction path. These 
conditions are very special, however. They serve 
to illustrate the constriction effect and the focus- 
ing effect free from the complications of a change 
in quantum number. In any actual surface, we 


in a circular channel. The vibra- 
tional component of momentum 
is changed on passing around the 
bend. 


should not expect to find the Schrédinger equa- 
tion separable in this fashion, and, even for very 
smooth reaction paths, the probability of a 
change in vibrational quantum number is quite 
high. This is demonstrated by the consideration 
of an energy surface consisting of three sections. 
A straight channel is joined to one of circular 
reaction path, which in turn is joined toa straight 
section. Consider a particle rebounding elastically 
from the walls of the channel. Let the last colli- 
sion in the straight section be just where the 
channel begins to bend. (See Fig. 2.) The total 
bend of the channel is 28+, and the last collision 
within the bend is y before the end of the bend. 
The path makes an angle a—y with the exit 
channel, whereas it made an angle y with the 
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Fic. 3. Rectangular corner. Classically, all particles 
whose trajectories lie below and parallel to PA are 
transmitted. 
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channel before the bend. This shows that the 
vibrational energy has changed in rounding the 
bend. 

Thus, in any actual case, there will be reflec- 
tion and transmission of the incident systems 
into states having a different vibrational quan- 
tum number from the initial one. The exact 
treatment of this problem is more difficult than 
when no shift in vibrational state is possible, be- 
cause the ease with which the shift occurs de- 
pends quite markedly on the peculiarities of the 
particular energy surface in question. Often the 
possibility of shifting into a new vibrational 
state makes the net reflection much lower than 
if this possibility were absent. In other cases, 
there is in reality very little reflection, but that 
which does occur is almost completely reflection 
into a different vibrational state. 

Some of the main features of these phenomena 
are illustrated by the behavior of systems travel- 
ing down a straight channel having a sharp right 
angled turn. The potential energy is zero within 
the channel and infinite outside it. (Fig. 3). 
Classically, the reaction is represented by a 
stream of particles rebounding from the walls of 
the channel. The component of momentum per- 
pendicular to the reaction path is proportional 
to the square root of the vibrational energy of 
the system. The component parallel to the reac- 
tion path is the translational momentum of the 
system. Hence, the trajectories of all systems 
having the same vibrational energy make the 
same angle with the reaction path. Let 


tan 0=),/p.=[Evin/(E—Evin)}! (13) 


where E= (p,?/2m) + (p,?/2m) is the total energy 
of the system. Then, referring to Fig. 3, consider 
a beam of particles crossing MN, each of which 
has the same total energy and the same distribu- 
tion of energy between vibration and translation, 
analogous to a single pure quantum vibrational 
state. The value of @ is then fixed by Eq. (13), 
but the position of the point P may be chosen at 
random so long as it lies on MN. However, in- 
spection of the geometry of the figure shows that 
all those particles whose trajectories are parallel 
to PA but lie above it will be reflected, while 
those whose trajectories are parallel to PA and 
lying below it will be transmitted. Calling the 
critical distance MP, d, and the width of the 
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channel /, we have, from the figure, 


l-—d l d 





2 sin , oe @ sing 
tan 6=3$+d/21l, 
derit =1(2 tan @—1). 


(14) 


The fraction of particles incident on MN whose 
trajectories are parallel to PA and which are 
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Fic. 4. Rectangular corner. Classically, all particles 
whose trajectories are parallel to and above PB are 
transmitted. 


transmitted is: 
T,=2 tan 0-1. (15) 


From Eq. (14) itisapparent that when 6=tan™ 3, 


no particles are transmitted. For smaller values 
of 6, inspection of the figure shows that there is 
also no transmission. When @=tan~ 1, all the 
particles are transmitted and for larger values of 
6, this is still the case. However, though the total 
energy and its distribution in vibrational and 
translational parts determines the value of tan 6, 
it does not determine its sign. We must now con- 
sider the transmission of particles whose trajec- 
tories are perpendicular to the line PA in Fig. 3. 
Figure 4 demonstrates the critical conditions for 
this case. All those particles whose trajectories 
are parallel to PB and lie above it are trans- 
mitted ; those lying below are reflected. The frac- 

tion transmitted is given by 
Funte pS ad ned (16) 

l l 
Again, for tan 0=3, T,=1. The over-all fraction 
transmitted will be the average of JT; and 7: since 
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the probability that 6 is positive is equal to that 
that it is negative. For the various ranges of 6 
we have the following results: 


T=3(7:+T.2)=tan@ O<tané@<1 


(17) 
T=1 tan @>1. 


If all vibrational energies are equally probable, 
the average classical transmission coefficient is 
obtained by integrating Eq. (17) with respect to 
Eyi»/E, using Eq. (13). The result of this graph- 
ical integration is 7=0.785. For comparison with 
the quantum mechanical treatment, the classical 
transmission coefficient was calculated by Eq. 
(17) for the case where E is restricted to integral 
multiples of E,yi,. The total transmission coeffi- 
cient for an incident particle with one quantum 
of vibrational energy is shown in Fig. 5. The 
partial transmission coefficients for transmission 
into states having one and two quanta of vibra- 
tional energy were computed from Eq. (17) and 
added to give the plotted curve. As in the quan- 
tum case to be discussed, the slower moving inci- 
dent particles have a greater chance of trans- 
mission than those having a greater proportion 
of translational energy, but there are no periodic 
critical values for which transmission is zero. 
The quantum mechanical picture replaces the 
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(b) 


Fic. 5. Transmission coefficients for a rectangular 
corner, as a function of the total energy of the system. 
(a) Particles incident in the lowest vibrational state. (b) 
Particles incident in the first excited vibrational state. 


particles by a train of plane waves whose in- 
tensity at any point is proportional to the density ) 
of matter there. The boundary conditions of Fig. 
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3 cannot be satisfied by any single solution to 
the Schrédinger equation. However, in each of 
regions marked J, JJ, and JJI such solutions can 
be found. The actual state of the assembly* will 
be represented in each region by a superposition 
of the possible stationary states, each with the 
appropriate amplitude. In region J, it is repre- 
sented by a term corresponding to an incoming 
plane wave in a particular quantum state, to 
which must be added an infinite train of waves 
traveling in the opposite direction in every pos- 
sible quantum state, including those for which 
the momentum is imaginary. Such a function is 


tk tpi. 


V,! =sin —y ex | - 7] 
‘ Tio eT 


“4 ™j tp; 
+> Ri; sin —y exp Ka (18) 
j=l l h 

where / is the width of the channel; ;, the trans- 
lational momentum of a system in the kth vibra- 
tional state; and R,; is the arbitrary amplitude. 
Those terms in Eq. (18) for which 9; is imaginary 
correspond, not to traveling waves, but to waves 
with exponential damping, such as are usually 
found in diffraction effects. The R,; are deter- 
mined by the condition that the wave functions 
and their first derivatives must be equal at each 
boundary of every two regions. Thus, we have 
found it necessary to introduce reflection into 
vibrational states differing from the initial one 
in order to satisfy the boundary conditions. 

The mathematical details are carried through 
in the Appendix. Here we shall point out some 
general results of this type of calculation. The 
most striking result is the oscillation of the trans- 
mission coefficient as the energy of the particle 
increases. Whenever the translational energy of 
the incident system becomes equal to the vibra- 
tional energy of some permissible quantum state 
for the transmitted particle, we find a zero in the 
transmission coefficient. Thus, when the system 
is incident in the first quantum state, its vibra- 
tional energy is one unit. If the exit channel is the 
same width as the incident channel, the vibra- 
tional states possible in it have energies of 1, 4, 


*The term “assembly” here denotes the collection of 
representative points which correspond to the equilibrium 
distribution of the reacting atoms in all their possible 
configurations. 





9, ---m? units. Hence a system incident in the first 
vibrational state will be totally reflected when its 
total energy is 1, 2,4,5,8--- (i.e., for any energy 
E=j?+n*, where n and j are integers). Figure 5 
shows the transmission coefficients for a particle 
having one quantum of vibrational energy ini- 
tially. The classical transmission calculated in the 
manner previously described is also shown. 

Reference to Fig. 5 shows that slowly moving 
particles are most easily transmitted, i.e., those 
whose total energy is only slightly greater than 
one unit. However, whenever there is a possibility 
that the transmitted system has translational 
energy, the transmission coefficient is markedly 
increased. Thus, at total energy equal to 4.5 
units, the coefficient for transmission into the 
lowest quantum state (no change in vibrational 
energy) is very small, in agreement with our 
previous observations for this case. However, the 
chance that the transmitted system be in the 
second quantum state is very much greater, so 
that the total transmission coefficient is greater 
than for systems of slightly less energy, e.g., 3.5 
units. We conclude that the presence of the 
corner tends to enhance the probability for a 
change in vibrational quantum number. 

The case may be considered as analogous to the 
transmission of electromagnetic waves. Waves of 
very short wave-length (or fast particles) travel 
in straight lines with negligible diffraction effects 
at irregularities in the reflecting surface. They are 
highly reflected at a sharp bend. Waves of 
longer wave-lengths (slow particles), such as ra- 
dio waves, are easily diffracted and can be trans- 
mitted around sharp corners with much less 
difficulty. Wave motion of all kinds, electromag- 
netic, elastic, acoustic, exhibits these same general 
transmission properties. 

It would be of interest to know to what degree 
the rapidity of oscillation of the transmission co- 
efficient is affected by the asymmetry of the po- 
tential surface, the removal of the discontinuous 
edges, and the averaging due to thermal equi- 
librium. For the case of a right-angled channel, 
calculations in the appendix show that narrowing 
the incident channel increases the rapidity of os- 
cillation, whereas narrowing the exit channel has 
relatively little effect. This is evident from the 
fact that total reflection with no change in quan- 
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tum number occurs whenever the total energy of 
the system is given by 


“9 2 
Ff .% 

fa—+— 
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(19) 


where / is the width of the incident channel, w is 
the width of the exit channel, and j is the vibra- 
tional quantum number of the incident system. 
The number n takes all integral values. For fixed 
j, diminishing / gives more closely spaced values 
of e than does diminishing w. 

It is also possible to calculate the transmission 
coefficient for an energy surface like the preced- 
ing, but in which the channels meet at an angle 
6 given by /=w sin @ where/ and w are the widths 
of the entrance and exit channels, respectively. 
In this case, the transmission coefficient shows a 
fluctuation with @ similar to that calculated for 
changes in energy. Thus, we should expect op- 
timum angles, differing for different energies, for 
which the transmission would be greater than at 
other angles. 

Some physical processes may be described as 
taking place on a potential energy surface on 
which there is total reflection due to the absence 
of a pass through the high energy region. When 
sound is absorbed in gases, the net result is the 
transfer of the translational energy of the sound 
wave to the internal vibrations of the colliding 
molecules. Such a process may be pictured on an 
energy surface having a blind alley. (Fig 6.) 
The reflection coefficient then becomes a measure 
of the ease with which the interconversion of 
energy takes place. By considering a straight 
channel of the type described above, one end of 
which is closed, we find that there can be no 
interconversion of energy if the channel is per- 
fectly symmetrical. However, if it is L-shaped, or 
T-shaped with side arms of unequal length, there 
will be an interchange. Again we find resonant 
energies for which there is no interchange just as 
in the cases of total reflection from an open-end 
channel. 

It is evident that the quantitative details of the 
behavior of the transmission coefficient are very 
sensitive to minor variations in the potential 
energy surface. Until these surfaces are known 
in detail and with greater accuracy than that now 
available we can hope to calculate no exact values 
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Fic. 6. Closed channel giving total reflection. 


for the actual transmission coefficients. However, 
the results obtained in the simple examples dis- 
cussed above indicate the general nature of the 
problem. Reflection on the potential energy sur- 
face can arise from three sources: energy humps, 
constrictions, and bends in the reaction path. If 
the hump or constriction is sufficiently gradual 
and the reaction path is straight, we have seen 
that there will be. little tendency for the system 
to change its vibrational state in traversing the 
activated state. However, if the reaction path is 
curved, there is a strong probability that the sys- 
tem will undergo a change in vibrational state, 
whether it be transmitted or reflected. The mag- 
nitude of this probability is an oscillating func- 
tion of the total energy of the system and of the 
angle through which the reaction path turns. It is 
more probable that the reflected or transmitted 
system be moving slowly than that it be moving 
rapidly, since the curvature induces transitions 
to the highest available vibrational state. The 
actual chance of reflection likewise oscillates with 
the energy of the system and the angle through 
which the reaction path turns. Slow systems are 
more easily transmitted than fast systems, in con- 
trast to the case in which there is no change in 
quantum number, where the opposite is true. 
However, the more highly energetic systems are 
less easily transmitted than those whose total 
energy is low. The analogy with the optical prop- 
erties of radiation is very marked and carries over 
even into the classical treatment, which is analo- 
gous to geometrical, or ray, optics. 

Any observable experimental result due to 
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these effects can be expected only if it appears 
that the peaks in the transmission coefficient are 
not obscured by the statistical averaging due to 
the fact that in practice we do not deal with 
systems all having exactly the same energy and 
exactly the same vibrational energy, but with an 
equilibrium mixture. Since in any actual chem- 
ical reaction, the initial systems are distributed 
statistically with respect to their total energies, 
the experimental transmission coefficient will be 
an average of those calculated here. Since the 
energy fluctuates in a range kT about the mean, 
we can estimate the average transmission coeffi- 
cient by taking it to be 


EqtkT Eo+kT 
c= [ (Beer | f eF/RTGE (20) 
E E 


o-kT o-kT 


where E, is the mean energy, x, the exact trans- 
mission coefficient and xk the experimentally ob- 
servable coefficient. If «(£) oscillates rapidly as 
E is increased by kT, the average value x will be 
a smooth function of E and the maxima and 
minima would not be observed. In Fig. 5, the 
energy unit is h?/8mi/*. For ortho-para hydrogen 
conversion, /=0.2 10-8 cm. Taking m to be the 
mass of the hydrogen atom, the energy unit is 
11.5 kcal./mole, or at 300°K, 19.2 kT. For 
heavier systems and for wider channels, the en- 
ergy unit is smaller. The peaks appear to be 
from } to 1 energy unit in width, so that averag- 
ing over a width 2kT would not smooth them out 
in this case. Removing the discontinuity and in- 
troducing asymmetries into the energy surface 
will narrow the peaks, so that it would be unsafe 
to predict any noticeable experimental deviations 
from a transmission coefficient of unity except in 
the case of very light atoms, or of isotope separa- 
tions. For these latter reactions, the relative 
masses of the isotopes fix the separation of the 
two systems on the energy scale of Fig. 5, so that 
whereas when one has an energy such that there 
is little transmission, its isotopic counterpart 
with the same energy may automatically fall on 
a region of Fig. 5 where there is good transmis- 
sion. Moreover, the angle of the bend is a function 
of the masses of the atoms in the system, and will 
be different for different isotopes. Since the trans- 
mission coefficient is also a function of the angle, 
there is a second possibility that different isotopes 
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have different transmission coefficients. These 
effects will be larger for light atoms, and in- 
deed, for thesé the observed separation factors 
are much larger than for heavier isotopes. How- 
ever, any quantitative prediction of these effects 
must await more exact knowledge of the energy 
surfaces involved. 

These results have several experimental impli- 
cations. In the first place, they indicate the im- 
possibility of measuring exactly the height of the 
energy surface by an experimental determination 
of the activation energy of a reaction. The virtual 
barriers due to curvature and constriction of the 
reaction path amount to at least 0.5 kcal. and 
prevent any closer estimate being made from ex- 
perimental data. Since all calculated energy sur- 
faces depend so far upon knowing the height at 
some point from experimental data, it is evident 
that they will not be useful in general for pre- 
dicting the second-order effects dealt with here. 
Although energy surfaces sufficiently accurate to 
enable the prediction of relative rates of reaction 
of isotopes are not known, enough factors have 
been shown to enter so that one might expect 
unusual separations in some cases. A third point 
of interest for experiment is that it seems likely 
that most newly-formed molecules exist in ex- 
cited vibrational states. This ought to be observ- 
able if a spectroscope were placed in the direction 
of flow of a streaming system of the reacting 
molecules. This would amount to observation of 
infra-red chemi-luminescence, and should give 
information on many hitherto unobserved higher 
vibrational levels. 


SUMMARY 


The quantum effects have been investigated 
for a number of idealized energy surfaces and 
have been found to be small except possibly for 
reactions involving isotope separation and for re- 
actions at low temperatures. It is shown to be 
highly improbable that no vibrational transitions 
occur during the course of a reaction, but that 
they must take place after the system has passed 
through the activated state. Consequently, the 
Eyring formulation of the activated state method 
is justified in its essential features. 
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APPENDIX 
Hyperbolic Channel 


Here we shall write down explicit expressions 
for the wave functions and the separation con- 
stant used in calculating the virtual potential 
barrier to which reflection may be ascribed. 

For our purposes it will be sufficiently accurate 
to take for the solution of the ® equation :4 


b= (1—k?* sin? ¢)-? 


2rR (?* : 
Xsin [a= f (1—k? sin? 6)'de+B | (21) 
RA do 
where 
Ak? [ k cos 2¢ 
4n?R%2(1—k? sin? $o)? 
sd i (22) 
16(1—k? sin? go)? 
k=(E/[A+E/2])!; A=h/(2mE)}. (23) 


Here A is the de Broglie wave-length of the inci- 
dent waves. This form for the solution of the 
Mathieu equation is similar to the asymptotic 
solutions of Langer and Jeffrey. It differs from 
their solutions in having the constant A different 
from unity. By using the value of A given above, 
it may be easily verified by direct substitution in 
the differential equation that ® satisfies equation 


Ph 82°mR? 
d¢? h? 


exactly at the point ¢=¢po. The error in the solu- 
tion may be estimated by the difference between 
the value of A calculated for the two extreme 
values of ¢ and that for the midpoint between 
the two hyperbolae. As the energy of the system 
becomes progressively larger, Eq. (21) becomes a 
better solution. For small values of the energy, 
we must confine our discussion to the small 
values of A. 

4 (a) R. E. Langer, Trans. Am. Math. Soc. 36, 637 (1934) ; 
(b) Harold Jeffreys, Proc. London Math. Soc. 23, 437 
(1923); (c) M. J. O. Strutt, Ergeb. d. Math. u. Grenzge- 
biete 1, No. 3 (Berlin, 1932); (d) Whittaker and Watson, 


Modern Analysis (Cambridge Press, 1927), fourth edition, 
Chapter XIX. 





A={1+ 








E 
(5 cos 26+n)@=0 (24) 
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The boundary conditions that ®=0 when 
= o— (4)A and when ¢= ¢+(3)A are satisfied if: 


27rR 





nyt = 


oth /2 
f (1—? sin? ¢)'d¢+B, 
0 


4 


IR oo-Al2 
ma =A—— | (1—k? sin? ¢)'do+B, (25) 
RA Jo 
IeR eootal2 
Nr = NT — Nox = A—— 
RA Joo-a/2 


x (1—R? sin? $)!d¢. 


For a given value of the energy, this equation 
determines the allowed values of k (and hence the 
separation constant, A). The integrals occurring 
in Eq. (25) are the incomplete elliptic integrals of 
the second kind and they have been tabulated 
many places in the literature.® It is therefore a 
relatively simple matter to compute the right- 
hand side of this equation as a function of k and 
to find those values for which the equation is true. 

If the angle between the asymptotes of the 
limiting hyperbolae is small, i.e., of the order of 
30° when the energy is not too large, it is possible 
to expand the integrals in powers of A: 


pot /2 
{ (1—k? sin? ¢)'dd=A(1 —? sin? ¢) 
o0—A/2 
A? / k? cos 2¢ k* sin? 2 
x{1-=( + — )|-.26 
24\1—k’ sin? 4(1—k? sin? dp)? 


Then it is easy to obtain the first two terms in the 
expansion of \ in powers of A: 

















po 50 24-9) (27) 
- —— cos 2¢9—A° 
8mR7A2 2. °° 
where 
nh? [1 Ecos2 EF? sin? 2¢ 
b(E) = fi = Po oe 0 


8mR°2A112 12n2E, 48n'E? 


PR I cos 260( ae ) 
2n?x* n° Eo 
——— ' $782 
x (1 yy sin? 60) +-(=) 


E 
X sin? 260(1 ote sin? +e) | (28) 
n*Eo 


Ey=h?/8mR?A? cos? ¢o. 

















(29) 


5 For example, see Jahnke und Emde, Funktionentafeln 


(Teubner, Leipzig, 1933), second edition, p. 124. 














Rectangular Corner 


Consider next a potential energy surface which 
consists of a rectangular trough of width /, 
within which the potential is zero, and outside of 
which it is infinite. We shall study the effect of a 
sharp right-angled corner in the trough. This 
potential surface may be divided into three re- 
gions, as shown in Fig. 3, region I being the inci- 
dent channel, region II, the square corner region, 
and region III, the exit channel. 

The Schrédinger equation is easily solved in 
regions I and II (Fig. 3). If we consider an inci- 
dent plane wave which represents the reacting 
system in its initial state in region I, superposing 
the reflected waves gives us the wave function: 
Wit =D (Rij exp [tpjx/h] 


J . 


+exp [ —ip;x/h ]6,;) sin 9. (30) 


In region III, there are only outgoing waves of 
the form 


YWM=> Ti; exp [ipjyy/h] sin <2. (31) 
j 


In region II there will be standing waves of the 
form 
. TT ; 
W.=>° Az; sin Pi sin a 
j h l , 


+ B,; sin “s sin oy, (32) 


The total energy of the incident wave is E= 
(p;?/2m) + (7?h?/8mil?). The first term is the trans- 
lational energy, and the second is the vibrational 
energy for the jth vibrational state. At each junc- 
tion of two regions, the corresponding wave func- 
tions and their derivatives must be equal. After 
making the expansions 


Me a 
sin —x=)) Cj, sin —x 
h 2 , l 
and 

. Pj . mn 

sin —y=)> Cyn sin —y 
h n l 

application of these boundary conditions gives 
us the following expressions: 


Rij= (4: sin H1—exp [—ipa/h Ji) 


Xexp [—ipd/h], (33) 
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PUR jexp [ipjl/h]—exp [—ipjl/h]éx,3) 
1 


vaiiaee cos Pinay (ile, (34) 
h hs i 
. Pi rr 
T).j;= By; sin Ys exp [—ip,l/h], (35) 
ip; . 
=, Tsexp Libil/h] 
1 
° . n 
= Bi cos 14 E (—1) Alay (36) 
Lt on 


These are to be solved simultaneously for R,;, 
Tx;, Axj, and B;;. The C,; are known functions. 
Putting Eq. (33) into Eq. (34), and Eq. (35) into 
Eq. (36), and rearranging we get: 


-" exp [—ipjl/h }(Axj+276:;) 
7 
TN 
=X (-1)"-BinCoj, (37) 
-" exp [ —ipjl/h |Bi; 
1 


== (-1) 


™m 
j ArnCnj (38) 





These equations must be solved for A;; and B,;. 
This is most easily done in matrix form. Writing 
the following matrices: 


b; rev mj 
re exp [—ipjl/h]éi;, Nejg=(— 1) bei 


Eqs. (37) and (38) become 


—(A+2i)P=BNC, —BP=ANC (39) 


whence 
A=-—21(1—G’)", (40) 
B=-—2i(1-—G*)"G (41) 
where 
G=NCP-, (42) 


Equations (40) and (41) in conjunction with 
Eqs. (33) and (35) complete the formal solution 
of the problem. 

The matrix elements G;.; are given by: 


Gij=(—1)**122°R7 sin gi 
Xexp [—1ge]/qi(qi?—m7?) (43) 


h 
eee qs=pd/h=a(8mPE/W— 7). (44) 
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The necessary matrix elements of (1—G?)~! had 
to be calculated numerically. Although this is 
an infinite matrix, the finite square 6X6 matrix 
was found to yield sufficiently accurate numerical 
values. 


Effect of Asymmetry 


To determine the effect of asymmetry of the 
potential energy surface on the transmission co- 
efficient, the previous calculation was repeated 
for the case where the width of one branch of the 
potential energy surface is a fraction, 1/a, times 
the width of the other. The wave function in 
region IT, Fig. 3, is then 


. vis ; 
WU=>> Ax; sin Pi, sin ~y 
j h l 


: ' 
T. 


+ B,.; sin aoe sin Pi, (45) 
w h 


where p;’ is now the momentum of a particle in 
the jth vibrational state in the exit branch of the 
channel. 

The analysis proceeds as before, and results in 
these expressions for Ry; and 7;,;: 


Rij= (Ax; sin ag;—exp [ —1aq; ]6x3) 
(46) 
(47) 


Xexp [—iaq;], 
Tx;= By; sin (g;'/a) exp [—7(q;'/a) ], 


where A;; and B;; are the respective elements of 
the matrices: 


A=-—2i(1-G"G’)—, 
B=—2i(1—G"G’)"G", 
2(—1)'tikj sin ax(8?—k?)} 

Xexp [ix(a?s?— 7*)*/a] 

(w/a) (a6? — 2) (a6? — ak? — 7?) 
2(—1)*+*k ja? sin [(a/a) (a?6?—k?)*] 
Xexp [ima(— 7*)*] 
m (a? — a?j”) 1(a26? — k? — a?7?) 


(48) 
(49) 


, (50) 








Gis= , (51) 
where &=8m/[E/h? and a=w/I. 

The general form of the transmission coeffi- 
cient curve is much the same as for the sym- 
metrical case, but the zeros in 7;,; are differently 
spaced. 
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For a=2, i.e., if the exit branch is twice the 
width of the entrance branch, the zeros are twice 
as dense as for the case a= 1. Higher values for a 
lead to still denser zeros in the transmission co- 
efficient. This means that as the total energy of 
the system increases the quantum mechanical 
transmission coefficient will oscillate the more 
rapidly the wider the exit channel is compared 
to the entrance channel. Since no actual system 
is composed of particles all having exactly the 
same energy, the average transmission coefficient 
will not display the maxima obtained for a single 
particle, unless the peaks are considerably broader 
than the energy range of the particles. This is the 
case for a=2, but for higher values the peaks 
become much sharper. 

Taking a=1 corresponds to passage from the 
wide channel into the narrow one. The general 
formulae remain the same, but the zeros are now 
spaced much as before. For a=} they come at 
E=4, 5, 8, 9, 13. Since a vibrational quantum in 
the narrow channel is four times the size of that 
in the wide one, there can be no transmission 
below E=4 when passing from wide to narrow 
channel. The subsequent zeros are the same as 
for the case a=1. 


Total Reflection 


The simplest energy surface for total reflection 
is a closed straight channel. In this case, however, 
no exchange of vibrational and translational 
energy can occur. This is immediately evident 
when the wave functions for incident and re- 
flected waves are equated at the closed end, x=0: 


Tj awk 
> Rij sin —y=sin —y. (52) 
j l l 
Hence, Ri.;=6.;. This is apparently the case 
whenever the energy surface is symmetrical about 
the reaction path. For an unsymmetrical surface, 
there will be interchange of energy. 

Such a potential surface may be idealized by 
considering a rectangular channel with a right 
angle bend in it, but closed at one end (Fig. 6). 
The general procedure is as before. In each of the 
three regions indicated in Fig. 6, solutions are 
found to the Schrédinger equation. In regions I 
and II, these are given by Eq. (30) and (45), re- 
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spectively. In region III we can have only 
standing waves, represented by 


VW Pt=> D,; sin ty 36) sin —s (53) 
j WwW Ww 


where g;'=;'w/h and p;’ is the momentum in 
region III. Again, the solutions and their first 
derivatives must be equal at the boundaries of 
regions I and II and regions II and III. Applica- 
tion of these conditions gives the following results 
for the reflected amplitudes: 


Rij= (4 xj Sin exp C—igaw/t us) 
Xexp [ —ig;w/1], 
Axj= { —21(1—G’G)“},,;, 
2kj sin (qx'l/w)[cos (q;'w/l) 
+i sin (g;w/l) | 
(q;w/l) [8° — (PR? /w*) — 7?) 
2kj sin (g.w/l) sin (q;1/w) 





Giy=(— 1)" 





Gis (—1)*+4 
nj=(—1) I 620? — (kw? /22) — 72] 


Xsin g;’(/—b)/w 
Qu=1( Fl? —k*)}; 
gx’ = 1 (Fw? — k*)}. 


&=8mE/h?; 
(58) 


When b=1, Eq. (57) becomes indeterminate. The 
analysis is then simpler, however, and we find 
that for b=], 


Rij= — 643 exp [—12q;w/1]. (59) 


Hence there is no interchange of translational 
and vibrational energies if b=/. Inspection of Eq. 
(54) shows that, as before, the reflection coeffi- 
cient is an oscillating function of the total energy, 
having zeros at the same points as for transmis- 
sion. The form of the matrix coefficient A;,; is 
different from that of the preceding cases, in that 
G and G’ are no longer symmetrical in a and 1/a 
as before. This however, cannot affect the general 
shape of the curve for R, but only the magnitudes 
of its maxima, which appear much as in Fig. 5. 


Obtuse Angle 


The effect of varying the angle of the bend in 
the potential energy surface can be calculated for 
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Fic. 7. Obtuse-angled bend. The corner region has the 
shape of a right triangle. 


obtuse angles using the wave functions: 


Wi! = Do (Rij exp Ligsx/l]+exp [ —igx/1 6x; ) 


_ mj 
Xsin—v, (60) 
l 


x j 


x ; 
Yu=> ( @xjsin + Anjcos *) sin 7? 
' ] 


rh 
sin ’ 
W COS > 





+ Bi; sin (61) 


w cos d 


~~ 
y= os Tx; exp (ig; §/w ] sin I (62) 
5 WwW 


J 


where 9j=pj/h, qj =pj'w cos ¢/h, qj" =pj'w/h, 
p;is the momentum of a particle in the jth vibra- 
tional state in the incident channel, whose width 
is 1; p;’, is the momentum in the exit channel of 
width w cos ¢; and p;”’, is the momentum associ- 
ated with a channel of width w. Region II is a 
right-triangular region. (Fig. 7.) The cartesian 
coordinates in region I are x and y; in region III 
they are &, 7, related to x and y by the formulas: 


£=y cos ¢—X sin ¢, 


n=y sin ¢+<x cos ¢. (63) 


The general method of procedure is now the same 
as in the previous cases, since wave functions and 
their derivatives must match at the boundary of 
each two regions. However, we have an addi- 
tional relation between the coefficients @;; and 
Ax; which arises from the invariance of the 
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Schrédinger equation to a rotation of axes about 
the origin. Thus, when W;!! is transformed into é, 
n coordinates, it must satisfy the boundary con- 
dition, W;11(¢, 0) =0. This gives the condition 


E Ga; sin (q;¢/w) sin (wjé/w tan ¢) 
"2X Aay con (:h/m) sin (we fG/eo ton ¢) 
— Buy sin (q/'t/w) sin (nj¢/w cot ¢) (64) 


connecting @,; with Ax; and B,;. 
The equations to be solved for Rxj, T:.;, Axj, 
and B,; are: 


Ri; exp [ig; cot ¢]+exp [—ig; cot $d; 


= A,; cos (q; cot ¢)+@; sin (q;cot ¢), (65) 


4 " 
[Re exp [ig;cot ¢]—exp [—ig;cot $ ]éx;] 
= 7 [ae cos (g; cot 4) —Any in (g;cot 4) ] 


(—1)"xn 
n wcosd 

X [@e; sin (qjn/w tan ¢) 
* +Axj 00s (qyn/w tan 6) sin (xjn/w) 


BinCnj, (66) 


+ By; sin (q;'n/w cot ¢) sin (x jn/w) 
=) Ti; sin (xjn/w), (67) 
i 
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Tj ; 
> ——L 4; sin (gjn/w tan ¢) 
j wtang 


+A;z; cos (qjn/w tan ¢) ] cos (7jn/w) 


U 


qi , : ‘ 
+ Bi;— cos (g;'n/w cot ¢) sin (xjn/w) 
W 
qi 
——[ Qi; cos (qjn/w tan ¢) 
Ww 
—A,,; sin (qjn/w tan ¢) J sin (xjn/w) 
Tj 
w cot d 





By ;sin (qj'n/w cot ) cos (7 jn/w) 
” 


-_ 3 Py sin (rjn/w) (68) 


in which the expansion 


sin (q;'y/w cos ¢)= >> Cyn sin (rny/I) 


has been made. Similar expansions may be made 
in the last two equations so as to express both 
sides of each equation as a Fourier series in the 
same interval. The resulting coefficients may be 
set equal term by term, yielding matrix equations 
for the unknowns as in the previous cases. Since 
the expansions are complicated, though straight- 
forward, no general information can be gleaned 
from the equations without extensive numerical 
computation. However, it is apparent that the 
transmission coefficient will fluctuate much as 
before with changes in energy and in the angle ¢. 
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A theory of the dissolution of three-dimensional cross linked polymers is presented. The 
existence of a definite melting temperature is discussed in terms of molecular constants char- 


acterizing the gel. 


HEN a certain critical amount of cross 
linking occurs in a polymer composed of 
long chain molecules, continuously linked three- 
dimensional structures known as gels occur.! 


1 Pp. J. Flory, J. Am. Chem. Soc. 63, 3083 (1941). 


One of the ways of distinguishing these three- 
dimensional structures from linear chain poly- 
mers is that at ordinary temperatures the gels 
will exhibit a limited swelling in substances that 
will dissolve linear chains of the same chemical 
constitution. 
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The existence of a definite temperature above 
which the swollen gel will be completely dis- 
rupted by the swelling pressure of the solvent 
will be discussed here. The case is especially 
clear cut when the cross linking bonds are 
weaker than the primary valence bonds along 
the chain, because in this case the disruption 
will occur at the cross linking bonds. In general, 
however, disruption will occur at the weakest 
recurring set of primary bonds in the gel struc- 
ture, and the cause of the periodic weakness in 
the gel structure might be due to absorbed 
impurities. For simplicity we shall treat the case 
where the cross linking bonds are the weakest, 
but the same considerations apply with slight 
modification to the other cases. 

The disruption of gels at a given temperature 
bears strong similarities to the melting of a 
solid. The solid-liquid transition occurs when 
the unfavorable heat of the transition is over- 
come by the entropy of randomness gained on 
melting, in other words AH=T,,AS. In the case 
of the swollen gel, below a certain critical 
temperature, there will exist two phases, a 
swollen gel phase and a solvent phase. The 
solvent will swell the gel a definite amount until 
the osmotic forces driving the solvent into the 
gel will be balanced by the elastic reaction of the 
gel. The elastic reaction of the gel is an entropy 
elasticity, the driving force being the tendency 
of the long chain molecules to spend more time 
in the numerous compact configurations than in 
less numerous extended ones. At a certain 
temperature, the unfavorable heat for breaking 
cross bonds will be counterbalanced by the 
increase in entropy resulting from the fact that 
the released chains can assume the more coiled 
configurations. 

Most melting processes occur very rapidly, 
in fact the melting rate is determined by the 
diffusion of heat. The rate of dissolution of gels 
at their melting temperatures is quite slow 
because of the large entropy and energy change 
involved in the process. The rate of dissolution 
is given by (k7/h) exp (—AF'/RT), where 
AF' is the free energy required to break the 
bond and, for the case of primary cross links, 
AF' is high. 

Gelatin gels, which exhibit a reversible sol-gel 


transformation, appear to be a good example of 
the considerations presented here. They are 
presumably composed of long chain molecules 
with a molecular weight of about 30,000. The 
cross linking bonds are probably hydrogen bonds, 
and the weakness of these bonds explains the 
low melting temperature. The gels are highly 
elastic and when stretched exhibit double re- 
fraction, and give x-ray patterns as does stretched 
rubber. Below their melting temperatures the 
gels exhibit a limited swelling in water. The 
melting is reversible, and by carrying out the 
operations slowly, the melting and the setting 
temperatures can be made identical. In fact 
the slowness of the dissolution is the chief 
reason why the similarity to the solid-liquid 
transition has been obscured. For this same 
reason it is hard to evaluate the experimental 
evidence obtained to date to decide whether the 
melting occurs continuously or discontinuously 
as the temperature is raised, but there certainly 
does exist some evidence which favors the view 
that an abrupt change is involved.* 

The possibility exists that the breakdown may 
occur over a small temperature interval, and 
that extraction of soluble materials and absorp- 
tion of heat occurs before the large scale breaking 
of all the cross linking bonds has been completed. 
More experimental evidence is needed on this 
point, but the author believes that for a fairly 
homogeneous gel the dissolution is a cooperative 
phenomenon involving the simultaneous break- 
ing of many cross bonds. The process is co- 
operative because every cross bond that breaks 
results in a gel with a looser structure capable 
of absorbing more solvent, and this more swollen 
structure will be more readily soluble than the 
original structure as will be seen in the forth- 
coming mathematical discussion. Cooperative 
processes of similar type give rise to sharp 
Curie points. The sharp melting points in gelatin 
gels appear to favor this possibility. 

The gain in entropy that will ensue when the 


_cross links are destroyed can be calculated on 


*For discussion of this and the other experimental 
findings on gelatin gels outlined here, and also references to 
the original literature, see Chapter XV of Colloids by E. S. 
Hedges, published by Edward Arnold and Company (1931), 
and Chapters 19 and 20 of Colloid Chemistry by R. J. 
oa published by Houghton Mifflin Company 

1939). 











292 


the basis of the statistics of long chain molecules. 
It is known that the probability that the distance 
between the two ends of a long chain hydro- 
carbon molecule lies between r and r+dr is 
given by?~4 
4 
P(r)r’dr =—— exp (—?°/d*)r’dr, (1) 
m3 

where X is the most probable value of 7. \ can 


be shown for a hydrocarbon of x links to be the 
following : 


A= (2//3)lov/n (2) 


(the most probable value is equal to ./% times 
the root mean square value). 

For hydrocarbons /)=1.54A. For rubber (cis 
structure) Wall> has recently calculated a value 
of 2.01,/nA, and for gutta percha (trans struc- 
ture) a value of 2.90\/n for the r.m.s. value. 

We wish to calculate the entropy change in 
the following process: the swollen gel breaks its 
cross linking bonds and the long molecules 
dissolve in solvent that they had imbibed. There 
is no appreciable change in the entropy of 
mixing of solvent and long chain molecules, the 
greatest contribution arises from the fact that 
the long chain molecules can now coil up. A 
contribution of R cal./degree per mole of chains 
(the communal entropy)® arises because the 
centers of gravity of the molecules are now no 
longer localized in the network structure but 
can move through the entire phase. 

In order to calculate the entropy of dissolution 
of the gel, it is necessary to make explicit 
assumptions about the structure of the gel. The 
most complete treatment of gel structure! indi- 
cates that the possible structures may be quite 
complex. For purposes of calculation and com- 
parison we shall consider two cases. In the first 
case the highly swollen network will be con- 
sidered to have the aspect of a three-dimensional 
fish net, and the possible interpenetration of 
network units will be neglected. It is to be 


2H. Eyring, Phys. Rev. 39, 746 (1932). 

3 W. Kuhn, Kolloid Zeits. 68, 2 (1934) ; 76, 258 (1936). 

4 A. Tobolsky, R. E. Powell, and H. Eyring, Chapter on 
“Elastic viscous properties of matter’ in Frontiers of 
Chemistry (Interscience Publishers, in press). 

5F. T. Wall, J. Chem. Phys. 11, 67 (1943). 

6 J. Hirschfelder, D. Stevenson, and H. Eyring, J. Chem. 
Phys. 5, 896 (1937). 
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realized that much looser structures are possible 
and undoubtedly are important when the cross 
linking is small. 

If in the swollen state of the polymer the 
average distance between network points is 7, 
then using the Boltzmann relation S=kIn P 
the change of entropy in going from the swollen 
state characterized by r, to the dissolved state 
where the distance between network units is X, 
the gain in entropy per chain of the network is 


(3) 


If communal entropy be added, then the third 
term vanishes. 

Let v be the number of cross linking bonds 
effective in the network structure of the swollen 
gel per gram of polymer. Then, if y is the 
number of chains out of each cross bond, the 
number of chains s per gram of polymer is yv/2. 
In many cases y=4 so that the number of 
chains is twice the number of bonds, that is, 
s=2p. 

If the distance between neighboring bonds in 
the swollen gel is 7, the volume can be seen to 
be 8vr?/34/3, inasmuch as the cross linking bonds 
can for simplicity be arranged on a diamond lat- 
tice. Therefore, for the case of y=4, r?= V/1.54», 
where V is the swollen volume per gram of 
polymer. If M is the molecular weight of a chain, 
m the molecular weight divided by the number 
of atoms in a chain, and N is Avogadro’s number, 
then \?= (4/3)1,?(M/m) and s= N/M. 

The gain in entropy per mole of ‘‘chains”’ is, 
therefore, 


AS =k[(r?/d?) —In (r?/d?) — 1]. 


3m(2V) 3m(2V)! 
AS= ( —In ). (4) 
4(1.54N1,3)§Mi = 4(1.54.N1,°)§M3 








In terms of the density p of the polymer and 
the volume fraction of the polymers, and if the 
swelling is sufficiently large to neglect the 
logarithmic term 


(5) 





| 3m(2/¢)! 
141.54.) EM pt | 
The heat of melting for the case of y=4 is 


equal to E/2, where E£ is the energy per mole of 
cross bonds. Hence the temperature of melting 








we \v 


\v 


of 








is given by 





3m(2/¢)! | (6) 


T.=E / 2R| 
4(1.54.N19*)§ Mp! 


Another possible treatment of the entropy of 
distortion of a network is that which has been 
recently used by several authors for the de- 
velopment of the kinetic theory of elasticity.’ 
In connection with the theory of limited swelling, 
Flory and Rehner*® have developed an equation 
for the change in elastic entropy on going from 
undistorted polymer to the swollen state. This 
is based on the assumption that in the swelling 
process the relative configuration of the network 
points are not changed, and the swelling merely 
magnifies the distance between network points. 
On this basis r/X appearing in Eq. (3) can be 
taken to be the cube root of the ratio of the 
volumes of the swollen and unswollen states 
(which is 1/¢ in our notation). Their expression 
for elastic entropy per mole of chains using a 


-tetrahedral model for the network is 


ASatastic = — $R[(1/¢)?'—1]. (7) 


By a slight extension of Wall’s® analysis for 
entropy of stretching one can derive the following 
expression for the elastic entropy per mole of 
chains: 


ASeiastic= —$RL(1/o)'—In (1/p)'—1].  (7’) 


The entropy for the dissolved state of the 
network per mole of chains neglecting entropy 
of mixing will differ from the undistorted polymer 
by Rin (i/¢) plus the communal entropy R. 
If ¢ is small, In 1/¢ can be neglected in com- 
parison with 1/¢. The entropy change on melting 


7 E. Guth and H. James, Ind. Eng. Chem. 33, 624 (1941). 

8P. J. Flory and J. Rehner, New York Academy of 
Sciences, Preprinted manuscripts on conference on High 
Polymers (Jan. 1943). A more complete paper on this sub- 
ject by these authors is forthcoming. The author wishes to 
thank Drs. Flory and Rehner for useful personal discussion. 
°F. T. Wall, J. Chem. Phys. 10, 131 (1942). 
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(i.e., going from swollen state to dissolved state) 
is, therefore [from Eq. (7) or (7’) ] 


AS=3R(1/9)!. (8) 


For the tetrahedral model y=4 so that 
AH=E/2 and 


E 


Tm = SRO ” 


On the basis of their theory for ASelastic, 
Flory and Rehner have derived an expression 
for ¢ in terms of the molecular weight M, the 
molar volume V;, of the solvent, and the density 
of the polymer. When this is substituted in the 
expression for 7, one gets 


E 
Tn = ' 
3R(M/2pVi)2" 





(10) 


This method of treating the entropy of dis- 
solution is based on the assumption that the 
swollen gel structure is essentially the same 
(except for a homogeneous dilation) as that of 
the undistorted polymer, and that the tetra- 
hedral cells interpenetrate to the same extent 
in both cases. It must be remembered, however, 
that the undistorted volume of the polymer is 
determined mainly by the van der Waals volume 
of the atomic groups (i.e., by energy considera- 
tions) and not by statistical considerations, as is 
proved by the fact that the density of polymers 
of the same chemical constitution is independent 
of molecular weight. This fact, as well as the 
possibility that interpenetration of the elements 
of the network becomes much less pronounced in 
the swollen state, limits the accuracy of this 
method of calculating entropy change. 

Experimental work on a series of gels with 
similar values of E, m, and J, with Eqs. (6), (9), 
and (10) used for comparative purposes, should 
yield useful information on the structure of gels 
in the swollen and unswollen states. 
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By making use of Polya’s theorem and the groups of syvmmétries of the geometrical models 
representing unsubstituted parent compounds, it is possible to include the various kinds of 
isomerism in a single discussion. The number of isomers calculated by this procedure is the 


total number of all types. 





SING group theory, Polya! solved a general 
problem in combinatory analysis applied to 
geometrical symmetry, and his results have been 
applied? to the calculation of the number of 
structural isomers in various kinds of organic 
compounds. However, the applicability of Polya’s 
theorem to isomerism is considerably more 
general than this. It can be used, for example, 
to demonstrate mathematically the unity be- 
tween the three types of isomerism (optical, 
geometrical, and structural*) usually distin- 
guished. Thus, employing the proper three- 
dimensional models, one obtains by this pro- 
cedure the total number of isomers of all kinds. 
We may begin by restating Polya’s theorem 
in an immediately applicable form. Let 3% be the 
group of symmetries (of order h) of the geo- 
metrical model representing the unsubstituted 
parent compound (e.g., methane, ethylene, ben- 
zene, etc.), there being p positions in the parent 
compound at which substitution can occur. Let 
P be a permutation group isomorphic to 3. Let 
Hj, j2---ip be the number of permutations in P 
which can be written as the product of 7; cycles 
of order one, j2 cycles of order 2, ---, and j, 
cycles of order p. By placing substituents (hy- 
drogen included) 1, ¢2, «++, gp» at positions 1, 
2, --+, p, respectively, one obtains a configura- 
tion. Two configurations (¢1, g2, --:, gp) and 
(y1', go’, ***, Gp’) are said to be equivalent mod 
if there exists a permutation 


a a 
s=(. ) 
11 12 ve dy 


in P such that gi;= gi’, «++, gip=@p’. Let Akiko: 


1G. Polya, Comptes rendus 201, 1167 (1935). 

2G. Polya, Helv. Chem. Acta 19, 22 (1936); Zeits. f. 
Krist. 93, 415 (1936). T. L. Hill, J. Phys. Chem. (in press). 

3 We do not include here structural isomerism of the type 
(CH;—O—CHs, C:H;OH). 
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be the number of non-equivalent (mod 3C) con- 
figurations when k; substituents of a first kind, 
k» substituents of a second kind, etc., are sub- 
stituted into the parent compound. Then 


> A kyko---X1*"'x9*2- oe 


=; Dd Hisis---ipf fot ++fpi?, (1) 
where 
fr=1teitiat- ++, fom=l+xi?+x2'+---, 
+, fp=Ut xP +--+. 


The summation on the right is taken over all 
permutations in P and the expression on the left 
is the resulting polynomial in x1, x2, ---. That is, 
the coefficient 


hittee 


of x1 "1x9 *2- oe 


gives the number of isomers possible when hk 
substituents of a first kind, ky of a second kind, 
etc., are introduced into the parent compound. 

As has been mentioned above, this theorem 
has been applied to the calculation of the number 
of structural isomers in certain cases. However, 
according to the above definition of non-equiv- 
alence, it is clear that two isomers of any kind— 
whether they be optical, geometrical, or struc- 
tural—correspond to two non-equivalent con- 
figurations and will, therefore, be counted if the 
proper model is employed. That is, two isomers 


. (of any kind) cannot be superimposed one on the 


other by means of rotations or reflections or both. 
Therefore Ax;k2--- is the total number of isomers 
of all types. In the examples previously worked 
out (in structural isomerism) the geometrical 
models chosen have been two-dimensional and 
hence optical isomerism has been excluded. If the 
geometrical model is the correct three-dimen- 
sional one, no such restriction is imposed. 
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This method may be used to “‘prove’’ alge- 
braically the existence of two isomers in the case 
of a single asymmetric unit or a single geometric 
unit. Rules for calculating the number of stereo- 
isomers in carbon chain compounds by a process 
of combining the asymmetric and geometric 
units have already been worked out.* However, 
when the optical isomerism arises in a ring com- 
pound and is the property, in a sense, of the 
entire ring rather than of particular asymmetric 
carbon atoms, the use of Polya’s theorem is suf- 
ficient in itself. 


OPTICAL ISOMERISM 


The Asymmetric Carbon Atom 


‘ : / 
In defining an asymmetric carbon atom alge- 


braically one might expect that the existence of 
only one element, the identity element, in its 
group of symmetries would be an adequate ob- 
servation. It is, however, only a necessary con- 
dition and not a sufficient one, since the group 
is derived from axes of symmetry, whereas planes 
of symmetry must also be considered. Thus, 
each of the tetrahedral structures represented in 
two dimensions by 


| i | 
—_-e b—C —a oi —a 

a i b 

(1) (2) (3) 


has more than one element in its group of sym- 
metries (each having at least one axis of sym- 
metry), while neither 


a 


| cae 

o——« nor i--~+ 
; | 

b c 

(4) (S) 


has an axis of symmetry or more than one element 
in its group of symmetries, although (4) has a 
plane of symmetry and (5) has not. 

Resorting to Polya’s theorem, the geometrical 
model to be used here is obviously the tetra- 
hedron. As is well known, a permutation group 


isomorphic to the group of symmetries of the 


*G. E. K. Branch and T. L. Hill, J. Org. Chem. 5, 86 
(1940); J. K. Senior, Ber. 60B, 73 (1927). 





tetrahedron has 12 elements: One permutation 
can be written as a product of 4 cycles of order 
one; 8 permutations can be written as a product 
of one cycle of order 3 and of one cycle of order 
one; and 3 permutations can be written as a 
product of 2 cycles of order 2. Therefore, Eq. (1) 
becomes 


Dd Akike-+-x1*xght- +» =e (fit+Bfs'fil+3f.?). (2) 


If we choose fy=1+%1+%2+%3, all possibilities 
are included since this allows the substitution of 
three different groups in addition to the hydrogen 
left unsubstituted. That is, the coefficient of 
X1X2x%3 in Eq. (2) gives the number of isomers 
when four different substituents (¢1, ¢2, ¢3, H) 
are bonded to the same carbon atom. Similarly, 
the coefficient of, say, x1%2, X2%3, X1Xex3", or 
x1?x2%3 gives the number of isomers for structures 
of type (4) above, etc. 

Putting f=1+x1+x%2+43, h=1+x%+x7+x7, 
and f;=1+2x+2*2+,' in Eq. (2), one finds 


Dd A kikoksxy*'xq*2x03*8 = 2x xox3 +H 1X 2x3" 
+ 1X 3X2? +%2%3%1? + etC., (3) 


all coefficients being unity except the coefficient 
of x1%2x3. This serves, then, as an algebraic con- 
firmation of the possibility of two isomers when 
four different substituents are attached to the 
same carbon atom, and of only one isomer in all 
other cases. 


Allene Type Compounds 


Optical isomerism of the type found in 


$2 G4 

should be mentioned here. The parallel (||) and 
perpendicular (1) symbols indicate the relative 
planes in which the two ends of the molecule lie. 
With this in mind, it is easy to show that the 
permutation group isomorphic to the group of 
symmetries of this structure (replacing all y’s by 
H’s) has 4 elements: One permutation can be 
written as a product of 4 cycles of order one; and 
3 permutations can be written as a product of 2 
cycles of order 2. Then, 


: A A kike++-X1"1x9*?- oe 


=i(fitt3fr’). (4) 
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Choosing again f;=1+%1+%2+4; and fe=1+x/ 
+%2?+x3’, 


D A kikoksxy"xo**x3"8 = 6x xorg +3x1x0x3" 
+3x 1x28 +x;!+etc. (5) 


The terms written on the right side of Eq. (5) are 
representative. The coefficients give the total 
number of isomers (structural and optical) for 
the type of substitution indicated. For example, 
if all substituents are the same (x,‘), there is one 
isomer. If all substituents are different (x ;x2x3), 
there are six isomers: 


a Cc a b 
, rai \ / 
C2z2C=C C=#2C=#C 
if ~ / *% 
b c d 
a b a d 
C2#C=C C=z2C2#C 
i \ a \ 
d Cc b Cc 
a d a Cc 
a / ™, rd 
C2=2C2°C C=C=#C 
Yd \ - * 
c b d b 


Ring Compounds 


Cyclopropane, cyclobutane, and cyclopentane 
may be considered planar and, therefore, the fol- 
lowing may be used as models: 


(1) (2) (3) 


where one set of substituents is undertood to be 
above the plane of the ring and the other set 
below this plane. The groups of symmetries of 
(1), (2), and (3) are easily obtained, and using 
Eq. (1), one finds, respectively, 


(a) DY Abske---xy*ixg*?- - - 
=43(f +3 fo3+2f5?), 

(b) S Abkske-+-xy*txq*2- + - i‘. 
=$(fi8+5fo'+2f.), 

(c) DY, Abjke---xy*txq*2--- 
= yo( fi? +5fe+4fs*). 
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As an example we may calculate the total 
number of isomers possible when » OH radicals 
are substituted into cyclobutane. Since all sub- 
stituents are to be the same, f,=1+4x, fo=1+2’, 
and fs=1+x*. Then, putting these in Eq. (6b), 
we have 
D> Apu’ = x8 +47 + 6x9 + 725+ 1304 

+7x?+6x?+x+1. °*(7) 


The coefficient of x" gives the desired result. 
Thus, the 13 isomers corresponding to the term 
13x‘ arise as follows: 


Structural: 


on on on on ow on ow on On 
L j l | i | i j rn j 
ion OW OW 
on o! o ou 
' 1 ' 1 ' 1 ' ' ' ' 
on 
4 + 3 + 1 + 1 = 13 


Optical: 4 + = 











There is some question about the geometrical 
arrangement in six-membered and higher rings. 
For any arrangement, of course, one can obtain 
the group of symmetries of the model and use 
the above method. For example, if a planar ring 
for cyclohexane is valid, 


; A kyko++-X1"'x9*2- rae 


_ Ta (fil? +7 fo8+2fs4+2 fe") 


as may be verified. 

Finally, it should be mentioned that applica- 
tions to rings containing one or more atoms other 
than carbon and to fused ring systems are also 
possible. In fact, as is clear from the generality 
of Polya’s theorem, any rigid geometrical model 
may be subjected to this type of analysis. 


(8) 


GEOMETRICAL ISOMERISM 


The group of symmetries of the planar con- 
figuration 


1 Y3 
: ee 
C=C 
eT 
$2 94 


is isomorphic to that of allene type compounds. 
Therefore Eqs. (4) and (5) are applicable (the 
total number of isomers being made up here of 
structural and geometrical contributions). For 
example, corresponding to the term 3x,x2x;° in 





1e 
of 
or 
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Eq. (5), we have the structures 


a Cc a b a Cc 


ee i ee le Me ale 
C=C C=C 
rw” f ™% / \ 


b Cc Cc Cc Cc 


ISOMERISM OF COORDINATION NUMBER n 


Up to this point we have discussed isomerism 
in connection with carbon compounds. The 
treatment.need not, however, be confined to 
carbon, nor, for that matter, to coordination 
number four. 

Structures are known in which coordination 
numbers from one to eight are exhibited, and 


Kimball® has listed the possible spatial arrange- 
ments. Using the group of symmetries of any 
such arrangement, Polya’s theorem leads to the 
symmetry formula [Eq. (1) ]. Thus the total 
number of isomers may be obtained for any sub- 
stitution in which the substituents are monova- 
lent coordinating groups. 

Marchi, Fernelius, and McReynolds® have re- 
cently discussed the stereochemistry of coordina- 
tion number eight, using models and the method 
of Lunn and Senior.’ 

5 G. E. Kimball, J. Chem. Phys. 8, 188 (1940). 

6 L. E. Marchi, W. C. Fernelius, and J. P. McReynolds, 
J. Am. Chem. Soc. 65, 329 (1943). 


7A. C. Lunn and J. K. Senior, J. Phys. Chem. 33, 1027 
(1929). 
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WO general developments of the kinetics of 
polymerization reactions have been re- 
cently given. One! involves the assumption of a 
steady state in regard to the concentration of 
active nuclei; the other? is more general. 
Copolymerization reactions have been treated 
in a few instances*® by introduction of certain 
ad hoc assumptions. In view of present interests 
it seems desirable to undertake a more complete 
analysis of these reactions. Because of the dif- 
ficulty of a general development, the analysis 
rests on the steady state method which is valid 
as long as the rate of initiation is negligibly small 
in comparison with that of growth and that of 
termination. Mainly termination involving a 
radical chain and a stable monomer is considered 
at present. A few of the results are outlined in the 
following. 


1E. F. G. Herrington and Alan Robertson, Trans. 
Faraday Soc. 38, 490 (1942). 

2R. Ginell and R. Simha, J. Am. Chem. Soc. 65, 706, 
715 (1943). 

3 R. G. W. Norrish and E. F. Brookman, Proc. Roy. Soc. 
A171, 147 (1939); F. T. Wall, J. Am. Chem. Soc. 63, 1846 
(1941). 


We postulate the mechanism: 


ee ky ky 
Initiation: Ny—— 1; No——nn1, 

ka1(r, 8) koa (r, 8) 
TPN ris; Nrs PN rs4ty (1) 








Propagation: 1,5 


— kgu(r, 8) kgo(r, 8) 
Termination: 2,,——>Ny41s3) Mrs— Neu 





N, and N2 are mole fractions of monomer 1 and 
2, respectively ; ”,, are mole fractions of polymer 
radical of specified composition. N,, are stable 


polymers. Since 
dn,,/dt=0, 


we have from (1) 
OrsNre=BriaMriastYrsiMrea 7,S=l, 
a9N10 = ky Ni, (1a) 
aio = ki2Ne, 

with the definitions 

Ors=[ka(r, S) +Rai(r, s) |i 


k q k: , No, 
Brs=kai(r, s) Mi, +[he(r, s)+Rs2(r, s) ]Ne2 


Yre=kan(r, s) No. 


The solution of the difference equations (1a) 
can be found exactly as a complex function of the 





298 H. 


two probabilities for growth of nuclei: 
, 
Wik=Bix/Oiny Wik =Vik/ ine 


This solution expresses the total number of ways 
a chain of composition (7, s) can be built. The 
probabilities may be expected to be functions of 
the composition of the unstable polymer. In 
view of the scarcity of experimental data, it 
seems hardly warranted to carry through an 
exact treatment. As a first approximation we 
may, therefore, replace the probabilities by 
average values obtained on substituting a mean 
composition for the individual chain composition. 
The error will be greatest in respect to the pure 
polymers, ”,o and mos. However, since their 
number is small, the error will not be large. With 
these assumptions the solution of (1a) reduces to 


ky r+s—1 ke r+s—1 
nn=arla')| *( )+z( ) @) 
ko Ss ko» r 


where k2; and ko are the average rates of propa- 
gation: 
ka Do nin= DL kali, k)nix (3) 
i,k i,k 
with a similar definition for kee. In this equation 
ko, and kee are independent of the composition 
of monomer in accordance with our assumption. 
The first term in (2) enumerates all ways of ob- 
taining the polymer 1,, from a nucleus 0; the 
second term, from a nucleus 79. 
Going to the kinetics, we have 


dN, , , 

= —kuNi—-M [Roi (t, k)+kai(t, k) |nix 
i,k 

and a similar equation for N2. Introducing aver- 

age rates of termination k3; and kz defined in 

analogy to (3), we find from (1a): 


kuNitkwN2 


P giaitattermete, 
* RNs + kaos 


i,k 
Substituting into the rate equations for N; and 
N2 we obtain for the over-all rate 


dz/dt= (koe —ka)2[ (Ru2z+hi2)/(Rsiz+hse) ], (4) 


where z= N,/N; and all rates of propagation are 
much larger than the rates of termination. This 
equation or its integral allows a determination of 
certain combinations of rate constants. It may 
be noted that the over-all course of the reaction 
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does not follow a simple first-order law as is 
found in simple polymerization.'? If the second 
species is used up more rapidly than the first, 
z—o for t>«; otherwise, 2-0. If zo, the 
limiting rate is 

1 dN, Ruka 1 dNe Rirkoe 

N, dt ka ' Ne dt ka 


which is a first-order law; similarly a limiting 
second-order over-all rate results for a second- 
order initiation. For the final number average 
molecular weight, we find with the above ap- 
proximations simply 


a M,N, (0) + M,N, (0) 
~ (es1/ ker) Ni (0) + (Fes2/ ee) No(0) 


M, and M, are the respective molecular weights 
of monomer; N,(0) and N,(0) represent initial 
mole fractions. 

In the case of radical-radical termination with 
a rate k3, the over-all rate can be shown to be: 


(Roe — ko1) 
(Rs)? 


For the change of size distribution with com- 
position we find with (1), (2), and (4): 
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From the definition of the w’s we see that they | 
can be written as functions of z. The equation | 
may then be integrated to give the final size dis- | 
tribution. Formulae for the weight average : 
molecular weight may also be derived in this | 
manner. 
These preliminary results are presented as a bis 
contribution to a treatment of copolymerization 
which will allow a quantitative evaluation of ( 
kinetic data and of the structure of copolymers. 
A detailed paper dealing with this problem and 
extensions is being prepared. 
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T is recalled that the van der Waals attraction 
between atoms may be represented in good 
approximation by a sum of interactions between 
pairs of atoms.'! This approximation is obtained 
by applying perturbation theory to the second 
order. The resulting interaction potential is of the 
order of Va/r® where V is the ionization energy, 
a is their polarizability, and ¢ is their distance. 

If the perturbation calculation is pursued to 
the third order, interactions between triplets of 
atoms appear. The order of magnitude of the 
interaction is 

Vai 


, 
112°%23°P31° 


where 712, f23, and 73; are the distances between 
the atoms. The magnitude and even the sign of 
the interaction potential depend upon the shape 
of the triangle with the sides 7j2, 723, and 73;. Our 
purpose is to investigate the character of this 
dependence. 

The term in the third-order interaction energy 
which depends on triplets of atoms is 





Ho ;H jw 
W,/"= : 
1 jx0(W;— W.)(Wi— Wo) 
1403 41 


where W;, Wi, and Ws are the unperturbed ener- 
gies of the intermediate states j and /, and the 
ground state, respectively. In the numerator Ho;, 
etc., are matrix elements of the perturbation. As 
the perturbing potential we introduce dipole- 
dipole interactions between pairs of atoms. We 
further make the assumption that the lowest 
states of the interacting atoms are s states. Then 
in the state designated by subscript j and also in 
the state designated by subscript / two atoms will 
be in p levels. One of the three atoms is in the 
same p level in states 7 and /. Each of the other 
atoms is in the lowest s level in one of the states / 
orj and in a p level in the other of the states / or 7. 





‘Pauling and Wilson, Introduction to Quantum Me- 
chanics (1935), Chap. XIV. 
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The dependence of the perturbation potential 
on the shape of the triangle of the three atoms is 
due to two reasons. Firstly, the three matrix 
elements [Ho;, Hj, and Hi contain the three 
factors, 1/ri2°, 1/re3°, and 1/r3;° (not necessarily 
in the order given). Secondly, the dipole-dipole 
interaction depends on the orientation of the 
interacting dipoles relative to each other and 
relative to the line connecting them. The orienta- 
tion of a dipole depends on which particular 
substate of the degenerate p level was chosen. 
Subsequently one must carry out the summation 
over all the substates corresponding to the vari- 
ous p levels that can occur in intermediate states. 

The summation can be carried out by a 
straightforward but lengthy calculation. One 
obtains 


we a COS Y1 COS Y2 COS y3+1 
> = ’ 


a 
112° 23°P31 





where y; is included by the sides 712 and rs; with 
similar definitions for yz and y3. C is a positive 
quantity of the order Va’. Its exact magnitude 
depends on the energies of p levels in the three 
atoms and on the magnitudes of the dipole 
matrix elements. 

It is of interest to apply the equation to a few 
special configurations of atoms. 

Case 1. The equilateral triangle. For y1=y2=73 
the interaction energy is 


11/8 
WW," =C 


1 19°f23°731° 


which is positive, and hence the interaction is 
one of repulsion. 
Case 2. A right triangle. For y:=90° the 
interaction energy is 
1 
W,”"=C 


112°%23°%31° 


and again the interaction is repulsive. 
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Case 3. Three atoms in a line. For y1=y2=0 
and y3;= 180° the interaction potential is 


—2 


WwW,” = C. 


712° 23°%31° 


which indicates attraction. 


Case 4. One atom far removed from the other two. 
For 131->1 102% 31123 and 73220 so that the inter- 
action energy becomes 





(—3 cos? yi +1) 
C ; 


31°F 19° 


Wo= 











